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We report a thorough theoretical study of the low temperature phase diagram of CS2CUCI4, a 
spatially anisotropic spin 5* = 1/2 triangular lattice antiferromagnet, in a magnetic field. Our results, 
obtained in a quasi-one-dimensional limit in which the system is regarded as a set of weakly coupled 
Heisenberg chains, are in excellent agreement with experiment. The analysis reveals some surprising 
physics. First, we find that, when the magnetic field is oriented within the triangular layer, spins are 
actually most strongly correlated within planes perpendicular to the triangular layers. This is despite 
the fact that the inter-layer exchange coupling in CS2CUCI4 is about an order of magnitude smaller 
than the weakest (diagonal) exchange in the triangular planes themselves. Second, the phase diagram 
in such orientations is exquisitely sensitive to tiny interactions, heretofore neglected, of order a few 
percent or less of the largest exchange couplings. These interactions, which we describe in detail, 
induce entirely new phases, and a novel commensurate-incommensurate transition, the signatures 
of which are identified in NMR experiments. We discuss the differences between the behavior 
of CS2CUCI4 and an ideal two-dimensional triangular model, and in particular the occurrence of 
magnetization plateaux in the latter. These and other related results are presented here along 
with a thorough exposition of the theoretical methods, and a discussion of broader experimental 
consequences to CS2CUCI4 and other materials. 



I. INTRODUCTION 

The spin- 1/2 nearest-neighbor Heisenberg antiferro- 
magnet on the two dimensional triangular (hexagonal) 
lattice is the simplest theoretical model for frustrated 
quantum magnetism. Without additional perturbations, 
the model is believed to order at zero temperature 
into a 3-sublattice coplanar ground state.'^'Hl However, 
one may expect a strong sensitivity to additional per- 
turbations to the isotropic triangular lattice Hamilto- 
nian. The material CS2CUCI4 provides an interesting 
example of a spatially anisotropic spin- 1/2 triangular 
antiferromagnet.^ For several years, neutron scattering 
magnetization, and specific heat measurement^^ on 
CS2CUCI4 have intrigued the community with unex- 
pected behavior. These experimental properties have 
been suggested by a variety of authors, including the ex- 
perimentalists themselves, to indicate exotic physics such 
as a spin-liquid ground state, unconventional "spinon" 
excitations at higher energies, and quantum criticality. 
Theoretical work on this material has been intense.^ 

An advantage of CS2CUCI4 is that the small exchange 
constants and high degree of magnetic isotropy allow for 
a fairly accurate determination of several of the largest 
Hamiltonian parameters, by comparison with the mea- 
sured single magnon dispersion relation above an approx- 
imately fully-polarized state. The approximate Hamilto- 
nian, including a magnetic field, determined by the ex- 
perimentalists in this way is 

ij i 

Here z, j are sites of a stack of triangular lattices (see 
Figs. fl]andl6]). The principle exchange interactions de- 



termined in Ref. fT6l are J = 0.374 meV on nearest- 
neighbor bonds parallel to the b axis, significantly smaller 
J' = 0.128 meV ~ 0.34 J on diagonal bonds in the b-c 
plane, and quite small J" = 0.045 J along vertical bonds 
between adjacent triangular layers. A Dzyaloshinskii- 
Moriya (DM) coupling was also measured along the di- 
agonal bonds, with Dij = —Dji = :^Da where the sign 
in specific directions are as indicated in Fig. [6| and 
D = 0.020meV= 0.05 J. 
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FIG. 1: Magnetic sites and exchange couplings in a 2D trian- 
gular layer in CS2CUCI4: on-chain bonds J (thick lines || 6), 
frustrating diagonal bonds J' (dotted lines). Stacked lay- 
ers are coupled by nearest-neighbor exchange J" . Crystallo- 
graphic a,6,c axes are indicated. 



The relative smallness of the couplings other than 
J suggests that one may perhaps fruitfully regard this 
Hamiltonian as one of Heisenberg spin chains along the b 
axis defined by J, which are weakly coupled together by 
the remaining interactions.^^ This point of view was 
validated in Refs. [T2HT5| where it was shown that much 
of the observed low energy ordering and the high energy 
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inelastic neutron scattering data on CS2CUCI4 can be cal- 
culated starting from known exact results for Id Heisen- 
berg chains. Indeed, numerical approaches in Refs IT9H25l 
showed that, due to frustration, the diagonal interaction 
may be as large as J' / J < 0.7 while still retaining ap- 
proximate quasi-one-dimensionality. 

While this approach has been quite successful for 
CS2CUCI4, most notably in directly confronting data 
without any adjustable parameters^ there remain some 
puzzling features in the experiments. One of the most 
striking ones is the drastic difference in the low tempera- 
ture phase diagrams of the material in magnetic fields 
aligned along the three different principle axes of the 
crystal. Though some aspects of these differences were 
explained in Ref. 14, based upon the "standard" model 
in Eq. ([T]), other glaring discrepancies remain. In this 
paper, we resolve these outstanding differences between 
theory and experiment by correcting the standard model 

of CS2CUCI4. 

It is important to emphasize that the corrections to 
Eq. ([1]) must be small, because the standard model does 
an excellent job in explaining a large volume of experi- 
mental data. The parameters in Eq. ([T]) were determined 
by high field measurements of single-magnon spectra,!^ 
which leave little room for doubt of their correctness with 
relatively small error bars. Moreover, the same model, 
used at zero and intermediate fields, is quite success- 
ful in reproducing the full inelastic neutron spectrum, 
containing b oth co ntinuum and magnon/triplon (sharp) 
contributions .l^^nH Nevertheless, in some field orienta- 
tions, entirely different low temperature phases are ob- 
served in experiment than are predicted by the standard 
model. Thus, we must somehow explain major qualita- 
tive differences in the ground states of CS2CUCI4 in a 
field by very small corrections to of no more than a 
few percent! 

A key message of this paper is that, indeed, the frus- 
tration and quasi-one-dimensionality of this problem can 
and do amplify tiny terms in the Hamiltonian to the 
point where they actually control the ground state. Sen- 
sitivity to small perturbations is of course an often-cited 
characteristic of frustrated systems. However, the ex- 
tent to which this sensitivity can be fully character- 
ized in the problem under consideration here is, to our 
mind, unprecedented. Using the methods of bosoniza- 
tion, the renormalization group, and chain mean-field 
theory (CMFT), we are able to distinctly identify the 
hierarchy of emergent low energy scales that control the 
very complex ordering behavior of the anisotropic tri- 
angular antiferromagnet, in a magnetic field and with a 
variety of very weak symmetry-breaking terms. 

This paper contains many results, and a thorough pre- 
sentation of the methods required to obtain them. To 
briefly summarize, we have determined the ground state 
phase diagrams for the ideal two-dimensional anisotropic 
triangular antiferromagnetic Heisenberg model, and for 
the model appropriate to CS2CUCI4, in all three distinct 
field orientations, over most of the range of applied mag- 



netic fields. In the former, we find spin density wave 
(SDW) and cone states, and in the SDW, a family of 
quantized magnetization plateaux. In the latter, we find 
several phases, including an incommensurate cone state, 
an commensurate coplanar antiferromagnetic state, and 
a second incommensurate phase, descended from the an- 
tiferromagnetic one. The occurence of these phases de- 
pends crucially on the field orientation, and matches well 
with experiments on CS2CUCI4. The associated phase 
diagrams in the temperature-magnetic field plane are 
shown schematically in Fig. [2] Details of each phase and 
its properties can be found in the appropriate section of 
the main text. 



(a) hWa 



(b) hWb 



T 








Cone 


\ FP 



T 






2 








'ft 








i 




AF 













(c) hWc 



(d) Ideal 2d model 



T 


1 
















ft 




IC 


fAF^FP 


Q 









T 








^1 SDW 


1 g \ FP 







. (JT 



FIG. 2: Schematic phase diagrams in the temperature- 
magnetic field plane for fields along (a) the a axis, (b) the 
b axis, and (c) the c axis, (d) Schematic phase diagram for 
the ideal 2d {J-J') model. Here we use the abbreviations: 
FP = fully polarized state; AF = commensurate antiferro- 
magnetic state; IC = incommensurate state; CAF = coUinear 
antiferromagnetic state; and SDW = spin density wave state. 
The shaded areas denote the regions in which Dzyaloshinskii- 
Moriya and exchange non-trivially compete. For these re- 
gions, we do not have reliable theoretical predictions at this 
point. 



One noteworthy highlight is that, remarkably, when 
the magnetic field is in the h-c plane, the spin correla- 
tions impugn the popular interpretation of CS2CUCI4 as a 
two-dimensional "anisotropic triangular lattice" antifer- 
romagnet. In fact, in this very wide regime, in the ground 
state, the spins are more correlated in the a-h planes, per- 
pendicular to the triangular layers, than they are within 
those layers! Taking into account these correlations is 
crucial to obtaining the proper low temperature phase 
diagram. They lead to an enhanced sensitivity to some 
very weak second neighbor and effective "biquadratic" 
interactions, which are needed to stabilize the antiferro- 
magnetic and incommensurate states mentioned above. 

The remainder of the paper is organized as follows. In 
Sec. In) we present some necessary background, includ- 



FIG. 3: A single triangular layer of the cone state, illustrated 
for a field along the a axis. Circles with arrows indicate the 
sense of precession of the spins, as one moves along the x 
axis. This is most easily seen intuitively by comparing every 
other spin, which compensates for the natural staggering due 
to the underlying Neel correlations of the Id chains. Note 
that in the cone state, all spins precess in the same sense on 
all chains within a plane. For fields along a, however, this 
sense alternates between successive vertical layers, owing to 
the staggering of D. Within a single layer, cone states for 
fields along other axes are identical to this one after a global 
spin rotation. 




FIG. 4: A single layer of the AF state, illustrated for a field 
along the b axis. The spins lie in a plane spanned by the b 
axis and a second axis within the a — c plane but at a non- 
zero angle to both the a and c axes. The component of the 
spin normal to b is antiparallel on successive even (or odd) 
chains, so that the patten has period Ay = 4 along the c axis. 
This corresponds to a doubled cystallographic unit cell in this 
direction. 



FIG. 5: A single layer of the IC state, illustrated for a field 
along the c axis. Circles with arrows indicate the sense of 
precession of the spins, as one moves along the x axis. This 
is most easily seen intuitively by comparing every other spin, 
which compensates for the natural staggering due to the un- 
derlying Neel correlations of the Id chains. Note that in the 
IC state, by contrast to the cone state, spins on alternate 
chains in the plane precess in the opposite sense. 



ble in this orientation, and we argue instead that the 
inter-plane coupling J" becomes dominant. A subtle hi- 
erarchy of energy scales (see Fig. |9| is exposed, which 
leads to the establishment of competing antiferromag- 
netic and cone phases in this case. Then, in Sec. [VTl we 
consider the field along the final principal axis, c, 



where 

an additional symmetry- allowed DM coupling plays a key 
role. It leads to the new incommensurate phase and an 
interesting commensurate-incommensurate phase transi- 
tion. Having established all the ground state phases, we 



discuss some experimental consequences in Sec. VII We 
give the explicit spin structures, describe the NMR line- 
shapes, which provide a telling confirmation of the the- 
oretical results, and establish the nature of the T > 
phase diagram. We conclude with some brief discussion 
in Sec. |VIII[ Several appendices present details of calcu- 
lations underlying some of the results in the main text. 



II. BACKGROUND 



ing the standard model Hamiltonian for CS2CUCI4, the 
low energy properties of a single Heisenberg chain, the 
results of a space group analysis of allowed DM interac- 
tions in CS2CUCI4, and a summary of the general scal- 
ing and Chain Mean Field Theory (CMFT) approach to 
studying competing interactions. In Sec. |III[ we deter- 
mine the ground state behavior of the "ideal" model of an 
isolated spatially anisotropic triangular Heisenberg anti- 
ferromagnet, often presumed to apply to CS2CUCI4. It 
clearly disagrees with experimental results for all three 
field orientations, as already pointed out in Ref. 14, In 
Sec. |IV[ we successfully apply the standard model to the 
case with a magnetic field applied perpendicular to the 
triangular plane, the DM interaction D playing a crucial 
role in reconciling the behavior with experiment. Next, 
in Sec. [V| we study the case of magnetic field along the 
b (chain) axis. The D term becomes rapidly negligi- 



A. Explicit hamiltonian and coordinates 

For the bulk of this paper, we will adopt a simple or- 
thogonal coordinate system, with x along the crystallo- 
graphic b direction, parallel to the chains, y along the 
c axis, perpendicular to the chains within the triangu- 
lar plane, and z along the a axis, perpendicular to the 
triangular planes. Moreover, we will adopt a simplified 
geometry, which respects the topology and interactions 
between spins, but does not precisely reproduce the ac- 
tual locations of Cu^+ ions. A description of the actual 
ion locations is given in Appendix [B] In our simplified 
geometry, the spins form a set of regular triangular lat- 
tices stacked uniformly along the z axis, with spacing 
1 between spins in the x and z directions, and spacing 
Ay = 1 between chains in the triangular plane. 

In this representation, the lattice Hamiltonian for the 
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standard model is 

^sm = ^JSx,y,z' ^x-\-l,y,z—h- Sx,y,z^ ^Hi-\-H2^ (2) 



xyz 



with 



Hi - ^^J'Sx^y^z • 

xyz 



^x+^,y+l,z) (3) 



and 



H2 — ^ ^ Sx^y^z ' ^x,y,z-\-l 



(4) 



xyz 



B. 



Review of low energy properties of Heisenberg 
chains 



We give a brief synopsis of known results on the low 
energy theory of the single Heisenberg chain in a field, 

Hid = J^S{x) ■ S{x ^1) - h^S'ix). (5) 



Here x are taken as integers, and we have taken the z axis 
in spin space along the field. In Eq. ([5| the magnetiza- 
tion M = j^S^{x) = {S^{x)) is conserved (A^ is the 
number of spins), and so it is convenient to work at fixed 
M. The equilibrium relation between magnetization and 
field M(h) is known from the Bethe ansatz solution, see 
Ref.i26i 

For any M less than full saturation, \M\ < 1/2, the low 
energy theory can be described in Abelian bosonization 
by a single massless free scalar field 6 and its "dual" (j) 
(related to the canonical momentum conjugate to 6) 



[0{x),(l){x')] = -ie{x-x' 



(6) 



with 9(x) the Heaviside step function. The low energy 
Hamiltonian is then 



Ho 



dx^^id.^f + id^ef), (7) 



where the velocity v depends on M, see Fig. 9 of Ref. [27l 
At a given M, the fiuct nations of the "longitudinal" spin 
component along the field axis are gapless at wave vec- 
tors kx = O^TT ± 2^, with 6 = ttM. Similarly, the "trans- 
verse" spin fiuctuations perpendicular to the field axis 
are gapless at kx = ±2(5, tt. The lattice spin operator is 
decomposed thereby according to 



(8) 



X) ^ e 



-i25x 



Here the scaling operators Sq^S^j^25 describe longitu- 
dinal spin fiuctuations, and 8^251 transverse ones. 
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M = 


M 1/2 
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TABLE I: Scaling dimensions of scaling fields associated with 
spin fluctuations in the one dimensional Heisenberg chain at 
magnetization M. The third and fourth columns give the 
scaling dimensions in the limit of zero and full polarization, 
respectively. 



These should be assumed to vary slowly with x (and 
time). Note that the operators 5^ (/i = z^±) do not 
mean the Fourier components of S^{x). They can be 
expressed in terms of bosonic fields as follows 



SS{x) 



2 

i 



-A,e 



-27ri0//3 



i/36> ±i27r0//3 



(9) 
(10) 

(11) 
(12) 



The parameter (3 is obtained by solving the integral 
equations, ^^ ^^ and the amplitudes ^1,^2 5^3 have been 
determined numerically in Ref. 31, We note that the 
above effective field theory describes the long-distance 
correlations of the spin chain, and is a good approxima- 
tion beyond some cut-off length ao- For "generic" values 
of the magnetization, this is of the order of lattice spac- 
ing; however, it diverges near saturation (|M| -> 1/2), 
where it scales roughly as the distance between fiipped 
spins (antiparallel to the field), ao ~ (1/2 — |M|)~^. 

The parameter /3 = 27rR is related to the "compacti- 
fication radius" R. At zero magnetization M = h = 0, 
the SU{2) invariant Heisenberg chain has 27tR^ = 1. In 
the field /3 and R decrease (see Fig. [12] of Appendix [A]) 
toward the limit 27rR'^ = 1/2 as |M| ^ 1/2. The scaling 
dimensions for general M are given in terms of R, and 
in the limits of zero and full polarization, are listed in 
Tabled! 

In addition to the scaling fields above, which appear 
in the expansion of spin operators, we will also make use 
of the spin current, for the component of spin along the 
field axis. This has the form 

r{x) = ^[S^{x)S-{x^l)-S-{x)S^{x^l)] 

-J^{x). (13) 

Like the spin density S§, this has a simple bosonization 
formula, 



J' 



(3J 



(14) 



We note that at M = 0, v j J — i\ {2 takes a simple value, 
but at M > the coefficient decreases continuously. 
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It may also be useful to connect with the limit of zero 
field M = ^ = 0, which is probably more familiar. Here 
the Hamiltonian has SU(2) symmetry. In this limit, the 
three operators Re[<S^], Im[5^], and Re[5|_25] become 
unified into the three components of the Neel field N 
(scaling dimension 1/2), while the other three operators. 
Re [5^], Im[5^], and Sq, become the uniform magneti- 
zation operator L = Jn + J/,, where Jr/l are the chiral 
spin currents (scaling dimension 1). The remaining oper- 
ator, Im[<S^_25] becomes the staggered dimerization field 
s (scaling dimension 1/2). 




C. DM-ology 

Here we present the DM terms which correct the stan- 
dard model, as allowed by the space group symmetry 
of the lattice. They are derived in Appendix |B] Since 
CS2CUCI4 is an S = 1/2 system with a single unpaired 
electron in a non-degenerate orbital, we expect that spin- 
orbit effects to be perturbative. In this limit, the leading 
effect is to generate DM terms on the same bonds on 
which exchange interactions are present, and which are 
proportional to both the exchange coupling on that bond, 
and to the strength of spin orbit interactions. As a con- 
sequence, we need to consider DM terms only on bonds 
with reasonably strong exchange, that is, the intra-layer 
triangular lattice bonds. These come in two types: the 
on-chain bonds and the diagonals. We denote the DM 
vectors on the former bonds by D and on the latter by 
D'^ in analogy to J and J' exchange couplings on the 
same bonds. 

The space group symmetry of the lattice determines 
the pattern of relative signs of the DM vectors on each 
of these bonds (see Appendix [B|) . We find the following 
form 

Hd = y ^ |-C^?y,2; ' ^x,y,z ^ ^x-\-l,y,z 
xyz 



'^■^y,z ' ^x,y,z X ^x-\-l/2,y-\-l,z 
'^■^y,z ' ^x,y,z X ^x-l/2,y-\-l,z^ 1 



(15) 



where 



Da{-lYa + D,{-lYc, (16) 
±D',{-lYa + D',{-lY+^b ± Di(-l)^+^c, 

(17) 



and D'^ = D \s the DM term in the standard model. The 
relative signs of the DM vectors are graphically shown in 

Fig.m 

One sees that symmetry allows five distinct DM cou- 
plings: D = D'^, D'^, D'^, Da and D^. Of these, only D'^ 
and Dc will be invoked in the body of this paper. The 
remaining three can be safely neglected, as explained in 
Appendix [Fj 




..M 3 




FIG. 6: Distribution of the DM vectors. Sites 1,2,3, and 
4 correspond to {x,y, z), {x -\- -\- 1, z), (x, y, z -\- 1), and 
(x + |, ^ + 1, 2; + 1), respectively, (a) The signs 0, refer to 
the direction of the DM vectors along a axis, (b) The arrows 
indicate the direction of the DM vectors along 6 or c axis. We 
use the convention (see Eq. ( 15 )) in which the first spin in the 



cross product Dij • Si x Sj is the one with smaller x for the 
on-chain bonds (i.e. it is on the left side of the bond), and 
it is the one with smaller y on the diagonal bonds (i.e. it is 
the lower of the two spins). Note that a different convention 
is used in Refs. 1161321 



D. Competing Interactions, Scaling, and Chain 
Mean Field Theory 



We assume that all inter-chain couplings are weak. In 
this case, a scaling analysis based on a perturbative RG 
treatment is appropriate. This proceeds in a standard 
way. One integrates out short-distance modes (i.e. small 
X or large k^)^ progressively reducing the large momen- 
tum cutoff from its "bare" value A (of order the inverse 
lattice spacing, which we take in turn to be 0(1)) to 
Ae~^, where £ G [0, 00] is the logarithmic RG scaling 
variable. Equivalent ly, the corresponding real space cut- 
off is L = ae^, where a is a microscopic length (which 
we take to be 0(1)). As we integrate out modes, the 
couplings themselves are renormalized. For each given 
coupling constant 7^, which appears in the Hamiltonian 
diS H = i^o + ' ■ --{-JdxjiOi^ we can define a dimensionless 
coupling = ji/{vAj), which is measured relative to the 
typical magnitude of the terms in the free bosonic field 
theory. Equivalently, division by both a factor of vA^ 
(the typical energy at this scale) and an additional fac- 



6 



tor of Ki (a typical inverse length) are needed to render 
7i, which is an energy density, dimensionless. To linear 
order, each dimensionless coupling "flows" according to 
the RG equation 



di^i = (2 - Ai)7i. 



(18) 



Note that the factor of 2 in this equation, which comes 
from the normalization by is equivalent to the 

space-time dimensionality of the (l+l)-dimensional field 
theory of the spin chains. In RG schemes in which space- 
time is rescaled to keep the cutoff fixed, this factor arises 
directly from that rescaling. We prefer to formulate the 
RG without rescaling in this paper, so that all lengths, 
times, energies, etc. are explicit. 

Many of the above operators may be relevant in the 
RG sense. This means that, with increasing ^, the di- 
mensionless coupling constants increase^ which implies 
A < 2 for that coupling. The RG is valid only so long 
as the largest of these dimensionless couplings remains 
small. Crudely, then, we may determine the length scale 
^ at which inter-chain coupling becomes significant by 
the point i = at which the first operator renormalizes 
to become of 0(1), where ^ = ae^ . The length ^ de- 
fines a correlation length, below which the dynamics is 
approximately one-dimensional. 

If this first "diverging" operator is unique, one can 
often identify the nature of the associated instability of 
the decoupled chains. In many cases, this can be done by 
dropping the other operators, and treating the remain- 
ing one by a type of mean field theory. Sometimes it 
can be treated in a semiclassical fashion. None of these 
approaches are rigorous, but they are eminently reason- 
able, and are likely to correctly predict the nature of the 
resulting state. 

From this reasoning, we see that both the bare mag- 
nitude and the scaling dimension (relevance) of the dif- 
ferent interactions are important in determining the low 
energy state of the system. To be more concrete, con- 
sider two candidate operators, Oi,(92. Their renormal- 
ized coupling constants obey 



7.(0 



V 



(19) 



Here, since A = A^^o is 0(1), we replace 7^(^ = 0) = 
ji{£ = 0)/v = ji/v. Setting 7i(^i) = O, an 0(1) con- 
stant, we obtain 



7i 



(20) 



We expect that operator 1 (2) is dominant if (^2) is 
the shorter length. Hence the boundary between the two 
regimes, in which one or the other operator dominates, 
occurs when ^1 = ^2, or 



7i 



2-A2 



(21) 



Note that, although the scaling of this boundary is de- 
termined by this argument, the precise location is not. 



due to the ambiguity of C. We can understand this con- 
clusion also from the scaling of the ground state energy 
density, Sq^ which obeys 



fo(7i,72)=5-%(^5^-^%^6^-^^). 
Choosing b = {^)~'^^^'^~^^\ we obtain 



(22) 



^0(71,72) = (^)' /(72/(7i/v) 



(2-A2)/(2-Ai) 



), 

(23) 

where /(A') is a universal scaling function. If there is 
a phase transition as 71 and 72 are varied, there must 
therefore be a non-analyticity at A* = A** in /(A'), for 
some A'*. This gives an equivalent condition to Eq. (21). 



To precisely determine the phase boundary, we need to 
know A'* (or O), which means we need knowledge of 
/(A'). Such scaling functions are generally determined 
by the full RG flows out of the scale-invariant theory, 
and not just perturbative data. 

To resolve this ambiguity, we turn to an approach 
which is equivalent to the former one at the scaling level, 
but which is more quantitative. A natural choice is the 
Chain Mean Field Theory (CMFT), in which inter-chain 
couplings are treated by a self-consistent Weiss-type de- 
coupling, using the exact solutions of perturbed but de- 
coupled individual chain problemJ^. Our CMFT ap- 
proach is described in detail in Appendix |Dj In principle, 
it can be employed to determine a full mean-field phase 
diagram. Here, we will mostly use it in more limited 
ways, as convenient. To address the ambiguity discussed 
above, we use the CMFT to compute a putative ordering 
temperature, Tj, for each channel driven by an operator 
Oi. The instability which sets in first upon lowering the 
temperature, i.e. with maximal Tj, is assumed to be dom- 
inant. Another application of CMFT will be to compute 
the magnitude of the ordering induced by a coupling 7^, 
at zero temperature. This will be useful in making quan- 
titative estimates of more subtle smaller energy scales, as 
we will see below. 



III. IDEAL 2D MODEL 

In this section, we consider the behavior of the ideal 
2d model, described by the standard model in Eqs. (TJS) 
with J" = D = 0. 



A. Continuum limit 

For this case, the technology of the previous section 
can be straightforwardly applied. We begin with the 
naive procedure of simply inserting the decompositions 
in Eq. (|8| into the microscopic inter-chain lattice Hamil- 
tonian in Eq. ([3| (We will assess the need to go beyond 
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this approximation later). Specifically, we have 



'y,z;7r-\-26 



(24) 



+e 



i26x 



It might appear that a slightly different formula should 
be applied for odd chains, since in those cases in our 
convention the x coordinates are half-integer rather than 
integer. However, the differences can be removed by con- 
stant shifts of and (f) for the odd chains, without any 



further changes. Hence we can uniformly apply Eq. (24) 
to all chains. 

Having dropped the D term, the behavior of the model 
at zero temperature is independent of the direction of 
the field, and is a function only of the magnetization M 
and the magnitude of J^ Inserting the decomposition 
of Eq. ([8| into Hi in Eq. ([3|, one need keep only terms 
which do not oscillate, the condition corresponding to 
momentum conservation. Using the slowly-varying na- 
ture of the scaling operators in x (but not in y and/or 
z)^ one may take a continuum limit in x by gradient ex- 
pansion to obtain the lowest non- vanishing terms of each 
type. One finds: 



Hi ^ J'Y^ldx [2M^ + 25,%^o5,\i,,;o + 2 sin J [S^z-.. -253^+1^,,^ +25 + h.c] 



[-^Sy^z^ndxSy^l^^.^ + h.c] + COS (5 



'^y,z;26'^y^l,z;25 + '^y,z;-26'^y^l,z;-26 + ^-^-^ | 



(25) 



Let us now assess the importance of each of the terms in Eq. (25). This is accomplished by ranking each of the 



terms in order of increasing scaling dimension, or equivalent ly, decreasing relevance in the RG sense. Formally, the 
most relevant term is the first (M^) one, which is a c- number constant and hence of dimension zero. Though it is 
a "trivial" constant (at fixed magnetization) and hence does not affect the dynamics of the system, it is indeed the 
dominant correction to the ground state energy of the weakly coupled chains. Being positive, it implies an increase of 
this energy with increasing M, and hence a suppression of the M{h) curve at fixed external field h. This is calculated 
in Ref ri4[ and reproduced in Appendix [Aj The result agrees very well with experimental data on CS2CUCI4. 

Consulting Table [l| one sees that of the remaining terms, the third and fourth terms are presumably most important. 
The third term involves S^^26 operators, whose scaling dimensions approach the minimal value of 1/2 at small 
magnetization. Because this term lacks any derivatives, it achieves nearly the smallest total scaling dimension (^ 
2 X 1/2 = 1) for small M. The fourth term contains a derivative (which adds 1 to its scaling dimension), but 
contains operators, whose scaling dimensions decrease from 1/2 towards 1/4 near saturation. Thus the total 
scaling dimension of the fourth term decreases from 2x1/2 + 1 = 2 at small M towards 2 x 1/4 + 1 = 1.5 near 
saturation. This makes it less relevant than the third term at small M, but more relevant than it near saturation. 
The remaining (second and fifth) terms have larger scaling dimensions for all values of the magnetization. 

We therefore drop these less relevant terms, as well as the constant contribution to the energy to obtain Hi H[^ 
with 



^1 — y^fdx |7sdw^^,^;7r-2(5^^+l,;g;7r+2(5 ~ ^coneSy^z^TT^x^v+l.z^TT + ^.C.j, 

y,z 



(26) 



with 7sdw = 2 J' sin ^, 7cone = J'f^. Equivalently, using 
the bosonization formulae in Eq. ([9|, one can rewrite H[ 
in sine-Gordon form, 



H[ = Mx|7sdwCos[27r(^^,^ -^^+i,^)//3] (27) 
y,z 

with 7sdw = J^Af sin J and 7cone = J^Alf3/2. 



B. Phases of (27) 



The names of these coupling constants have been cho- 
sen to reflect their probable consequences. For small 
magnetization, where 7sdw is most strongly relevant, one 
expects collinear "spin density wave" (SDW) ordering of 
spins along the z (field) axis, with 



(28) 



which minimizes the SDW interaction term. Here G 
[0, 27r] can be arbitrary for each since the layers are 
decoupled. Near saturation, where 7cone becomes more 
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relevant, one expects a "spiraling" order of the compo- 
nents of the spins transverse to the field, 



(29) 



with some qo > 0^ and where = ±1 and are in- 
dependent for each z. Some incommensurate pitch qo is 
preferred by the derivative in Eq. (26), but is expected 



to be small as it is disfavored by the single chain Hamil- 
tonian. Because of the non-zero magnetization along the 
z axis, the spins in this phase sweep out a "cone" as one 
moves along the b (chain) axis in real space. 



Following the logic in Sec. IIP the dominant interac- 
tion, at each magnetization, is the one whose putative 
ordering temperature is largest. The estimated ordering 
temperatures for the SDW and cone states, calculated 
from CMFT, are shown in Fig. [7[ One can see that the 
two curves cross at M ^ 0.24, which therefore separates 
a region of SDW state at lower magnetization from a cone 
state at higher magnetization. 



0.4 
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0.0 
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FIG. 7: (Color online) Putative ordering temperatures within 
chain mean field theory for SDW interaction, (Tgdw, dotted 
(red) line) and cone (Tcone, solid (blue) line) interactions, 
within the ideal two-dimensional model. Dashed blue line, 
emanating from M = 1/2 point, represents Tcone oc (1/2 — M) 
(Eq. (31)) which describes the cross-over between Id behav- 
ior obtained by bosonization with CMFT, and 2d one, which 
takes over for M > 0.45. The larger ordering temperature 
is expected to be physical, while the putative lower tempera- 
ture transition will be suppressed by the dominant competing 
order. 

This change from an SDW to a cone state is pri- 
marily due to the variation of scaling dimensions with 
magnetization. As M increases, the spin correlations 
are increasingly XY-like, i.e. enhanced transverse to 
the field and suppressed parallel to it. While the SDW 
term is obviously more relevant than the cone term near 
M = 0, due to the derivative in the latter, the change 
in scaling dimensions eventually compensates. Equat- 
ing the two scaling dimensions, one finds R = Rc, with 
27rRl = (V5 -l)/2^ 0.62 (the golden ratio!). This oc- 
curs at magnetization Mc ~ 0.32 (i.e about 65 percent of 



the saturation value). This approximately recovers the 
more accurate estimate Mc ~ 0.24, obtained above. Sim- 
ilar estimate holds for the SDW-to-cone phase transition 
in a single zig-zag ladder.!^ 



C. Critical temperature 

Further details of the behavior of the critical temper- 
atures in Fig. [T] can be understood physically. The SDW 
critical temperature increases from small values at small 
M, to a maximum around M = 0.1, above which it again 
decreases. The initial rise is due to the partial relieve 
of frustration of the inter-chain J' coupling by increas- 
ing incommensurability of the longitudinal spin correla- 
tions. The ultimate decay of T^dw is due simply to the 
suppression of longitudinal spin correlations as the chain 
becomes more XY-like. The same increasing XY-like be- 
havior leads to the growth of Tcone with M. 

The two endpoints, M ^ and M ^ 1/2, require spe- 
cial consideration. Approaching zero field, the dominant 
SDW interaction vanishes. This case requires a subtle 
analysis of fluctuation-generated interactions, which was 
undertaken in Ref. TT, There it was observed that correc- 
tions to the naiVe continuum limit are crucial to obtain 
the correct behavior, which is neither an SDW state nor 
a spiral, but rather a commensurate, collinear, antifer- 
romagnet (CAE). This CAE state replaces SDW as the 
ground state near the M = limit, see Fig. |2]i, which is 
the reason for Ts^w/cone curves in Figure [7| (as well as in 
most other Tc vs. M figures in the paper) start not at 
M = but at a finite M = 0.02 value. We will not go 
into further detail on this point here, but simply mention 
that another instance of fluctuation-generated couplings 
will be encountered later in Sees. [VllVIl 

The other limit, M 1/2, can be attacked differently. 
At M = 1/2, one has full spin saturation, and the state 
is unique and trivial. Single spin-flip magnon excitations 
can be found exactly including the effects of J^ One may 
obtain in this limit a cone state as a magnon condensate, 
as in Refs. I35|36l In this formulation, it is clear that 
the critical temperature for the ordering must vanish as 
M ^ 1/2. However, this is not observed in Fig. [?[ This 
is due to a non-commuting order of limits. In the vicinity 
of saturation, scaling (see Sec. D6) in fact predicts that 
all physical quantities are functions of the combination 



1 

2 



(30) 



when ^ — M <C 1. The bosonization analysis carried out 
above is valid for S <C 1. However, for S 1 a different 
behavior obtains. Specifically, the critical temperature is 
expected to scale (up to logarithmic corrections) accord- 
ing to 



J' /J 



M 



(31) 



9 



where - for S <Cl and - S for S > 1. One 
can see that this form indeed vanishes on approaching 
saturation. The maximum of Tcone should be obtained 
by differentiating, occurs when J-'(S) = SJ^'(S), which 
imphes S of 0(1). Hence the maximum Tcone occurs very 
close to saturation, where ^ — M ~ J' / J ^ and its pres- 
ence is not captured in the bosonization result plotted 
in Fig. [Tj Hence the temperature Tcone is overestimated, 
leading to an underestimate of the magnetization of the 
crossing point from Fig. [?[ This effect is probably small, 
however, since this occurs relatively far from saturation. 



D. SDW and magnetization plateaus 

The above considerations treat only the most relevant 
terms in the effective Hamiltonian. This, however, ne- 
glects some important physics in the SDW phase. In 
particular, it misses commensurability effects, when the 
SDW period can become "locked" (i.e. fixed over a re- 
gion of field and temperature) to a multiple of the lattice 
constant. 

Microscopically, this effect arises from "umklapp" pro- 
cesses, which distinguish quasimomentum from true mo- 
mentum, violating absolute conservation of the former. 
The X component of the quasimomentum is important 
here, and hence umklapp events carry momentum 27r. 
Since the SDW carries momentum tt ± 2(5, a umklapp 
event occurs when a number k of SDW quanta are ab- 
sorbed or emitted, adding to ±27r. This condition is rig- 
orously derived below, where we outline symmetry con- 
siderations which fix the form of the allowed microscopic 
umklapp Hamiltonian completely. 



1. Symmetry constraints 

We start by analyzing how (py{x) transforms under dis- 
crete lattice symmetries.!^ For that, we re- write (Isl) as 



-Ar sin[^^^M -{tt- 2S)x]. (32) 

It then follows that translation along the chain trans- 
forms 4>y{x) as 

c^y{x) ^ c^y{x + 1) - - 25), (33) 

ZTT 

while translation along the north-east diagonal {y ^ y -\- 
1, X ^ X + 1/2) changes it to 

^y (x) ^ (x + 1/2) - ^ (^ - 2S) . (34) 

47r 

Spatial inversion {x —x) changes it as well: 



In addition, (j)y{x) is defined modulo /3 so that 

(l)y{x) (j)y{x) + P 



(36) 



must be respected also. 

We now specify the general form for the /c-th order 
umklapp term: 

^umk = ^ j dx tk{y) C0S[^ (l)y{x) + UJk]- (37) 

y 

where tk ~ 0{J) is the bare amplitude and ujk is yet 
undetermined phase. The periodicity requirement, Eq. 
( [36| ), implies that k must be an integer. The translation 
in Eq. ([33| changes the argument of cosine in Eq. (37) 
into 27rk(j)y/f3 -\- — k{7r — 26) which implies that 

k{7r - 2S) = 2iiu. (38) 

Since 5 = ttM, the above equation implies that allowed 
values of the magnetization are given by 



2v 



(39) 



where v and k are positive integers. This condition is 
equivalent to the magnetization quantization condition 
for a single spin chain. However, we will see that v 
and k are not arbitrary in the two-dimensional triangular 
lattice. 

The remaining symmetries, translation along the di- 
agonal of the triangular lattice, Eq. ([34| , and spatial 
inversion, Eq. (35), require that tk{y) = (— l)^^t/c and 
uok = — 7r/2 for odd k and cj/c = for even k. As a result, 
the most general form of the umklapp term, consistent 
with lattice symmetries and involving single chains and 
no spatial derivatives, reads 



(fc=odd) 
umk 



(/c=even) 
umk 



J2 / dxtfe(-irsin[— 0,(x)],(4O) 

y 

^ J dxtk{-i)y^cos[—u^)].{ii) 



2. Energetic constraints 

To proceed, we bring out the energetics associated with 
the underlying SDW order by making the shift 
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(42) 



SO as to minimize the leading SDW term in Eq. ( |27[ ) 
(see Appendix Dl). In terms of the transformed fields, 
the SDW ground state corresponds to an x- and y- 
independent configuration of (l)y{x) = ^o- This shift 
modifies Eq. (pOl) as 



(/c=odd) 
umk 



<Pyix) -^'i:- 4>y{-x). 



(35) 



y 



dx t/c(-l)^^^+^^ 



Ilk 27ik 
sml^^l co^\^^(j)y{x)\. 



(43) 
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We immediately conclude that odd-k umklapp processes 
must have odd v in order to be able to gain some energy 
- otherwise the sum oscillates with y and does not give 
an extensive contribution. Similarly, Eq.(41) becomes 



(fe=even) 
umk 



(44) 

implying that now v must be even. 

This interesting result allows us to finally represent all 
allowed umklapp terms in a single compact equation 



^ ^ dx ikCos[^(l)y{x)], (45) 



where tk = tk{cos[^] + sin[^]), and Eq. (|45]) must be 
supplemented by the important constraint 

iy = k (mod 2), (46) 

i.e. ly must have the same parity as k. 

3. Allowed plateaux 

It is natural to consider the plateaux in order of in- 
creasing k - this will be directly related to the robust- 
ness of the plateau (see below). One finds that the first 
non-trivial possibility (different from the non-polarized 
M = or the fully polarized M = 1/2 limits), corre- 
sponds to = 1, /c = 3 when M = 1/6 = |(^). Assuming 
for concreteness i^^s > we find that = /5(2n + l)/6 
(with n = 0, 1,2) minimizes Eq. (|45l (and, by construc- 
tion, the SDW interaction in Eq. (|27|)). Working back- 
wards through the chain of the transformations we find 
that Eq. (32) predicts (recall that (tt — 2S) = 27r/3 here) 



{S'yix))M=i/6 = M + A^i-l)y cos 



2ttx 7r(2n + l) 



(47) 

This equation describes the famous^^ up-up-down (uud) 
spin configuration of the 1/3-magnetization plateau, i.e. 
with two-thirds of the sites having a larger spontaneous 
magnetic moment than the remaining third. It also cor- 
rectly predicts relative arrangement of down-spins on 
neighboring chains: the system gains energy by cou- 
pling every down-spin with a pair of up-spins on adja- 
cent chains. The resulting pattern has down-spins in the 
centers of hexagons formed by up-spins. For the other 
sign, ik=3 < 0, one finds instead of the uud configura- 
tion one in which two-thirds of the sites have a smaller 
spontaneous moment than the remaining third. This cor- 
responds to the "quantum" magnetization plateau sug- 
gested in Ref. 001 where the magnetic unit cell is com- 
posed of a spin singlet on a pair of sites accompanied by 
an up-pointing spin. 

Other possible plateaux include M = 3/10 = |(^) 
{k = = 1) and M = 5/14 = f(^) {k = 7,iy = 



1). Importantly, several of the smaller- /c plateaux are 
excluded due to the 'mismatch' between the parities of k 
and u numbers. These include /c = 4, z^ = 1 which leads 
to M = 1/4 (one-half plateau) and k = = 1 which 
would result in M = 1/3 (two-thirds plateau). 



4. Effective two-dimensional sine- Gordon model 

We now use the RG to derive an effective two- 
dimensional sine-Gordon model. Our starting point is 
given by the following Hamiltonian: 



plateau 



^ j dx{^{d:,(t)yf - cos[y(0^-^^ 

-^Ox(Py{x) + tk COS[-^(/)J|, 



(48) 



which incorporates the shift Eq. (42). Observe the ap- 



pearance of the new, linear in spatial derivative term, 
which is added here^^ to describe variation of the mag- 
netic field h near the optimal plateau value h^^'^\ The 
optimal field is defined by the condition that the magne- 
tization in the absence of the umklapp term, Mo(/i), is 
given by the plateau's value, Mo(/i^^'^)) = M^^'^). Then 
g(/i)=27rA:(Mo(/i)-M(^'^)). 

We then iteratively integrate out high-energy modes, 
reducing the momentum cutoff from initial Aq ~ 1 to 
Asdw = Aoe~^^^^, as described in the Appendix [e) The 
new, reduced cutoff Agdw is determined by the condition 
that the renormalized SDW coupling, 7sdw(Asdw/Ao)^^'^'^ 
becomes comparable to the contribution of the gradient 
term to the energy density at the same scale, '^A^^^. 
Here Agdw = 2/(47ri?^) is the scaling dimension of the 
SDW cosine in Eq. (|48|. This leads to the estimate 
Asdw/Ao - (5^)i/(2^sdw), At this scale, SDW cou- 
pling is of the order vA? ~ {itdwl^^'"^^)'^'^^'^^ 
SDW term should be minimized. Therefore, the argu- 
ment of cosine is small and we can approximate 



1 



m{x,y)f. (49) 



The umklapp term has a scaling dimension, A/^ = 
(27r/c//3)^/(47r) = k'^ /{AttR^), which grows quadratically 
with the /c, and is thereby strongly suppressed by high- 
energy fiuct nations: 



(-^sdw) 



4(0)(^)- 



(50) 



At this stage it is convenient to define a rescaled field, 
cp = 27Tk(j)/P + 7rO(f/e), which includes a shift to achieve 
a definite sign of the umklapp term {Q{x) is the Heavi- 
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side step function). This gives the two-dimensional sine- 
Gordon Hamiltonian of the k-th plateau, 



plateau 



vq{K)_, 



,ipy{x) - I t/e I COS [(/:)] I . (51) 



Here u = v{l3/2iikf, Cy - ^^-^(x^^^/v^-dw) 2-A,dw . 



5. Plateaux width 



At T = 0, the two-dimensional sine- Gordon Hamilto- 



nian Eq. (51) can be analyzed classically.^^ We review 



this standard analysis as it is important both here and 
in Sec. IVI CI The classical sine-Gordon model exhibits 
two phases: commensurate, describing the plateau, with 
{dx^) = 0; and incommensurate, with (dx^) 7^ 0. The 
incommensurate state, which describes the SDW phase 
with field-dependent ordering momentum, is achieved for 
sufficiently strong \q\ > Qc. The critical value Qc is de- 
termined by the condensation of kinks^ when £^kink = 0. 
Here the kink represents the solution of Eq. ( [51] ) inter- 
polating between two degenerate minima of cosine po- 
tential: (p{x — —00,7/) = and (p{x = +00,?/) = 27r, for 
all y. One immediately observes that the linear deriva- 
tive term in Eq. ( [sT] ) contributes —vq{h)/k to the kink's 
energy (per unit length in the y direction). The rest fol- 
lows from standard steps, which show that energy of 
the kink, relative to the energy of the uniform plateau 
state with, for example, (p = 0^ is given by 



kink 



8\u\tk\ 



v\q\_ 
k 



Thus 




(52) 



(53) 



Since in the relevant range of the magnetic field dMo / dh 
is constant, the plateau width in field units, ^/^(/c,^/) is 
directly proportional to Qc. 

Focusing on the 1/3-magnetization plateau we can es- 
timate, with the help of Figure 12] that 27rR^ 
M = 1/6. This leads to (5/^(3,1) 



(J7^)9/2 



3/4 at 
The next 



most robust plateau is at 3/5*^ of the saturation magne- 
tization {k = = 1), for which ^/i(5,i) ~ [J' jv)^^!^. 
The existence of this plateau is unclear, since this mag- 
netization is close to the boundary of the SDW state, and 
indeed, the calculation in Figure [7| predicts that it falls 
outside the stability range of the SDW state. It is, how- 
ever, inside the SDW phase as estimated from the pure 
scaling dimension criterion 27tR^ > (a/5 — l)/2. Thus this 
plateau still seems a reasonable candidate for observation 
in some anisotropic triangular materials. 

Other plateaux, such as, for example, the one-half 
magnetization one (M = 1/4), are much narrower due 



to the equal parity requirement^ Eq. ( |46[ ), which implies 
k = S^i^ = 2. Then (5/i(8,2) ^ {J' making it very 

hard to observe indeed. These arguments make it clear 
that the 1/3 plateau (M = 1/6) is drastically more robust 
than others, and thereby it is expected to be much more 
commonly observed. We also reiterate that it should 
persist all the way down to the decoupled chains limit, 
J' = 0. 

It is interesting to compare our findings with those in 
Refs. ,42 43 which studied magnetization plateaux in a 
single zig-zag ladder, made of two spin chains (with ex- 
change J) coupled in a triangular (zig-zag) fashion by 
exchange J' . This geometry can also be viewed as a 
single chain with first and second-neighbor interactions 
Ji = J' and J2 = J. In this one-dimensional system a ro- 
bust 1/3-plateau is found to exist in the intermediate ex- 
change region 0.487 < J / J' < 1.25. In particular, it does 
not seem to extend far into the J' ^ limit, although 
a very narrow sleeve of the plateau phase cannot be reli- 
ably excluded by the current numerical studies. We note 
however that the much reduced extent of the long-range 
ordered plateau region, in comparison with the quasi-2d 
predictions above, simply refiects the reduced stability of 
the crystalline (more specifically SDW in this case) or- 
der in Id systems (at T = 0). In more technical terms, 
the magnetization plateau requires pinning of both the 
'center-of-mass' and relative combinations of 0i^2 fields 
in a two-chain system which, in turn, requires significant 
modification of the chain Luttinger parameter K from 
its bare value of 1/2, in notations of Ref. 43 , by various 
marginal (density-density type) inter-chain terms. Such 
modifications generically require J' ~ 0(1), which is the 
reason for the absence of the plateau in the J' — )► limit 
in this case. 



6. Critical behavior of the wavevector 

Our description can be extended to the neighborhood 
of the plateau-SDW transition, where the ordering mo- 
mentum \q{h)\ > qc shows abrupt variation with mag- 
netic field. Near the commensurate state, the incommen- 
surate SDW phase can be understood as a solit on lattice^ 
with a finite linear density of solitons^'^ In the di- 
lute limit UsW <C 1, where w = ^/u/\tk\ is the width 
of the soliton, the solitons repel each other, with an ex- 
ponentially decaying potential /7e~^/^. The pre-factor 
U = 32^/u\tk\ can be obtained by calculating the energy 
of two solitons separated by a distance x with the help of 
(51). As a result, the energy density of the dilute soliton 
lattice is given by 



sol.lat. 



(54) 



Here the last term represents the repulsion between the 
nearest solitons of the lattice. The optimal concentra- 
tion n*, for q > qcj follows by minimizing Eq. (|54|) (with 
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logarithmic accuracy): 



w Inl 



u 1 

27Tv{q-qc) J 



(55) 



This imphes that the shift of the ordering momentum 
SQ = 27rn* from its commensurate value Qk,u = '^TTu/k 
inside the plateau vanishes with an infinite slope, accord- 
ing to 



SQ 



-27r 



^ln|(/i-/i(^'-))/(5/i(,,,)|' 



(56) 



where the last expression is written with logarithmic ac- 
curacy. 



with experiments on CS2CUCI4. At zero magnetic field, 
the ground state is actually an incommensurate spiral, 
similar to that predicted at high fields in the previous 
section, rather than the collinear state produced by the 
fiuctuation-generated interactions. Moreover, the zero 
field spiral ground state appears to continuously evolve 
on increasing fields along the a axis, with no intervening 
phase transition before reaching a fully polarized ferro- 
magnetic state at the saturation field. The SDW state 
predicted in the previous section is entirely absent. 



IV. STANDARD MODEL: FIELD (INCLUDING 
ZERO) ALONG a AXIS 

The predictions of the last section for the idealized 
model with J" = D = unfortunately do not agree 



This behavior, however, is readily explained by the 
standard model including the D and J^' terms, as dis- 
cussed in Ref. 14. To proceed, we again apply the de- 
compositions, Eq. (24), to the these terms. Following 



the logic of the previous section, we keep only the most 
relevant contributions. This gives 



-25) 



+7^-2 {^y,z;T^-25^y,z+l]TT-25^ 

with 7^' = = 2ii = J". 



-Ai5x 



' ^y,2;7r+2(5^y,2+l;7r+256'^*'^'^) + 7± (^^2;7r^y,2+l;7r + ^y ,z;T:^t,z+l]T^) }' (^'^) 



At zero field, we can simplify, using SU(2) symmetry 
and ^ = 0: 

H'^ = Y.jdx[- Di-ir (iV+.iV-+i,. + N-,N+^,^.) 

+ j"Ny^,-Ny^,+ ^y (58) 

A. Competition between D and J" 



In Eq. (58) it is evident that both D and J" induce 
strongly relevant perturbations, with scaling dimension 
1. These are more relevant than any terms naively 
present in zero field (i.e. the two terms studied in the 
previous section), and much larger than the fiuctuation- 
induced correction (~ {J'Y / J^) with the same scaling 
dimension, which drives the formation of the collinear 
antiferromagnetic statd^in their absence. They also be- 
come more relevant with increasing field. Hence we ex- 
pect that these terms should control the actual ordering 
in CS2CUCI4 for this field orientation. 

However, it is not so obvious which of the two is dom- 
inant. Indeed, they actually compete. This can be seen 



r 



as follows. The D term is minimized (for D > {)) hy 
configurations in which 



(59) 



where the classical phase i^iay depend upon z. For 
such configurations, however, the J" term oscillates in 
sign with and hence averages to zero. To instead min- 
imize the J" term, one requires configurations in which 



{Ny,.)j" = i-irNfly, 



(60) 



where the unit vector fiy = (n^, Uy^riy) may depend upon 
y. For all such configurations, the D term vanishes when 
averaged over z. 

The balance of this competition is determined by the 
relative magnitudes of D and We rely again on 

the chain mean field method, which indicates that the 
DM term dominates for the parameters of CS2CUCI4 - 
as shown by the fact that the associated mean-field or- 
dering temperature in Fig. [S] is largest. This is essentially 
due to the fact that the DM interaction acts on twice as 
many bonds as does the interlayer exchange. 

A physical distinction between the two candidate 
states, which may be compared to experiment, is in their 
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FIG. 8: (Color online) Same as Fig. [t] but with the chain 
mean field ordering temperature due to interlayer interac- 
tions (Til, dot-dashed (magenta) line) and DM interaction 
(Td, dashed (cyan) line) included. Note that the highest two 
temperatures are Td and Tii, with the former substantially 
larger. Thus for a magnetic field along the a axis, where the 
D term is present, it is expected to dominate the ordering ex- 
cept in the region close to saturation, M > 0.45 or so. That 
later region, which is described by a 2d scaling, is represented 
by the shaded rectangle. 



vector chirality, which is concentrated on the diagonal 
bonds of the triangular lattice. Define 



(61) 



This quantity is non-zero in both phases. In the contin- 
uum limit, one obtains 



1 



y,z,± 



y±l,z 



h.c.) 



(62) 



Let us compare this chirality for the states favored by D 
and J" . For the D term, one obtains Xy,z,± ^ (— 1)^A^^, 
which is constant in the triangular planes but alter- 
nates between layers. This staggering of chirality along 
the crystallographic a axis is observed experimentally 
in zero field. For the J" term, one obtains instead 
X± ~ ^^(^^^^+1 + ^^^^+1)7 which can vary within the 
triangular planes but is the same in every such layer. 
Thus experiment supports the T)-induced order but not 
the J" one, in agreement with the calculation described 
above. 

Now consider non-zero field. In this case, returning to 
Eq. (57) with J 7^ 0, we see that the third term is os- 



cillatory, and can be dropped. The second term is less 
relevant than the first and fourth, and thus is also sub- 
dominant. One is still left with a competition between 
the first and fourth terms, again controlled by the bal- 
ance of DM and inter-layer exchange. Since several terms 
that had formed part of the latter coupling for h = are 
now removed, we should expect that the DM will be rela- 
tively enhanced and continue to win the competition for 
all values of the field. This gives a natural explanation 



for the continuity of the ordered phase across the range 
of fields observed in experiments (with this field orienta- 
tion). It is also interesting to note that the scaling di- 
mension of the Sy^^.^ operator decreases with increasing 
field, actually making the DM coupling more relevant. 
This is reflected in an increasing critical temperature for 
the ordered phase with increasing field, up to a maximum 
which occurs relatively close to saturation. Again, con- 
sulting Fig. |8) we see that indeed the D term dominates 
the ordering except very close to saturation (M > 0.45 
or so). In that narrow field window, the larger bare value 
of the cone interaction {0{J')) is sufficient to overcome 
its larger scaling dimension. 



B. Incommensuration of the ordered state 

Naively, it would appear from the above analysis that, 
in the region M < 0.45, a commensurate ordered state 
is induced by the D term. However, while the D term is 
indeed dominant in this regime, we still need to take into 
account the subsidiary effects of the J' interaction. It 
turns out that the cone coupling 7cone does not actually 
compete with so that it introduces a weak incommen- 
suration in the ordered state. 



To see this, we apply the expectation value in Eq. (59) 



to the cone interaction in Eq. ( |26| ) (using S^^^.^ = TV^^ 
in zero field), assuming is a slowly- varying function of 
x: 



(63) 



The (average) phase is obviously the classical analog of 
Oz^ so that there is a gradient cost obtained from Eq. ([7|), 
which should be added to the above term to obtain 

y,z 

(64) 



This is easily minimized with respect to i?^: 



>7c 



'-j3N'^{-iy 



qo 
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(-ir. 



(65) 



Here qo is the induced incommensurability. For zero mag- 
netization, we have qo = {AJ' / J)N'^ . However, Eq. (65) 
also applied to M > 0, if is replaced by the magni- 
tude of the expectation value of S^^^.^. In both cases, the 
result is that 



(66) 



where 0^ is arbitary for each layer, since we have up to 
now neglected inter-layer coupling. Let us compare to 
what is expected in the cone-dominated regime, M > 
0.45. Here we should apply the ansatz in Eq. (29), which 



minimizes the cone interaction, to the D term in Eq. (57). 
One obtains 



y,z 



dx{-DM\-ira, 



(67) 
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which is minimized by taking = sgn(Z))(— 1)^. Sup- 
pose that D is positive, which is always possible if we 
redefine z. Then Eq. (29) becomes 



(5+,; J = (-l)^1^|e 



(68) 



Comparing Eq. (|66[) and Eq. (68), we see the forms are 



identical. Thus, the regimes M < 0.45 and M > 0.45 are 
actually smoothly connected, and distinguished only by 
which interaction controls the largest part of the ordering 
energy. This means that D and J' do not really compete. 
Indeed, in Appendix |D 4[ we show that in CMFT both 
interactions together increase the critical temperature of 
the cone state. 



C. Interlayer correlations 



The expression in Eq. (68) contains an undetermined 
phase, Qz^ for each layer. One may look to J" ^ which 



has been neglected in obtaining the form in Eq. (68), to 



fix these phases. However, at the naive level of first order 
perturbation theory, this is not the case. In particular, 
taking the expectation value of the J" term in Eq. (|58) 



or the corresponding j'^ term in Eq. (57), one finds a 



oscillating result, which vanishes upon summation over y. 
This indicates that the effects of J" on the undetermined 
phases is second order in 

Such second order effects can be considered as a 
fluctuation-induced interaction, which can be derived in 
a similar way as in the calculation of Appendix [Ej One 
obtains 



AH 



dx cos( 



(69) 



with ferromagnetic Jlj ~ {J"Y /v > 0. Taking its expec- 
tation value, this terms splits the large phase degeneracy, 
leaving only two undetermined values. 



6, 



Oo for z even 
Bi for z odd 



(70) 



Some bare microscopic second neighbor exchange might 
contribute to J2 , but experiments indicate that the net 
result remains ferromagnetic, as there is no enlargement 
of the unit cell in the a {z) direction. 

For the standard model, the two remaining phase de- 
generacies are protected by symmetry. The overall U{1) 
phase, 60 + 6i is of course expected to be arbitrary, 
owing to rotational symmetry of the Hamiltonian about 
the a (z) axis. The relative phase, Oq — Oi, is pro- 
tected by translation symmetry, x ^ x + 1, under which 



V. FIELD ALONG b AXIS 

In this section and the next, we will discuss the physics 
determining the ordered ground states when the mag- 



netic field is normal to the crystallographic a axis. These 
cases are much more complex than above, because, as we 
will see, the ordering is determined by several distinct in- 
teractions which are important at different energy scales. 
The "cascade" of energy scales, which must be considered 
in turn, from largest to smallest, is indicated graphically 
in Fig. [21 



Energy 




£2d 



£3d 



£CIT 



FIG. 9: Cascade of energy scales operative for fields in the 
b-c plane. Symbols adjacent to the arrows indicate the inter- 
actions responsible for the associated (partial) ordering. The 
cascade can be quantified by different "condensation" energy 
densities, which give the lowering of the energy density due to 
the establishment of the associated partial order. At the high- 
est energies, between J and £2d J^^IV^I^, the system exhibits 
one-dimensional fluct uati ons. Here IV^I is the amplitude of the 
order parameter, Eq. (75). Between e2d and esd, the spins or- 
der within a-b planes, but the planes are not registered. Below 
£3d ^ Max{J2|'0|^, ^"^J iV^l^ji full three-dimensional order de- 
velops, which may be of antiferromagnetic (AF) or cone type. 
In the former case, there may be yet another lower energy 
scale, £ciT (which is of the same order as Ssd), below which 
the magnetic structure is completely determined. This may 
be commensurate or incommensurate, the latter occuring only 
for fields along the c axis, and is driven by DM interactions. 



A. Irrelevance of D term 

Having understood that the DM interaction D dom- 
inated the physics for fields along the a axis, we first 
consider its role in this field orientation. Remarkably, 
the change in orientation has a drastic effect. With a 
field in the b-c plane, the D term always involves one 
spin component parallel to the field and one perpendic- 
ular to it. Consulting the decomposition of spin opera- 
tors in Eq. (24), we immediately see that the dominant 



fluctuations of these two spin components are always in- 
commensurate. As a consequence, in the continuum limit 
all terms directly arising from D oscillate with a e^^*^^ 
factor. This makes them formally strongly irrelevant. 
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More physically, upon coarse-graining over length scales 
shorter than 7r/(5, these terms average to zero. Thus, 
provided that D does not scale to strong coupling under 
the RG before this scale is reached, they become negli- 
gible. This will be true everywhere except the low field 
limit. Specifically, since it has scaling dimension 1, the 
renormalized D term at this scale is of order D{7t/S). 
Demanding this be small compared to we obtain the 
criterion 6 = ttM > J for it to be negligible. 

Thus for most of the magnetic field range, we can drop 
the D term. This gives a simple reason why the exper- 
imental behavior in CS2CUCI4 for this field orientation 
is completely different from that with field along the a 
axis. It is tempting to expect instead that the ideal 2d 
model discussed in Sec. [in| should apply. A comparison to 
experiments strongly suggests, however, that this is not 
the case, except perhaps at high fields. Unfortunately, 
no published low temperature NMR or neutron data are 
available in this field orientation. However, magnetiza- 
tion and specific heat measurement show a single phase 
occupying the region below the high field cone state and 
above the low field DM- dominated spiral. By compar- 
ison to the same measurements along the c axis, which 
are strikingly similar, it appears most likely that this in- 
termediate phase represents a commensurate state. In 
the ideal 2d model, we would be forced to interpret it 
instead as an incommensurate SDW. This seems unten- 
able, as such an SDW phase should also show a 1/3 mag- 
netization plateau (see Sec. HID), of which there is no 
sign. Thus we conclude that the SDW phase has been 
superseded by another competing state. In the follow- 
ing, we find an explanation for this competing state as a 
commensurate "antiferromagnetic" (AF) phase. 



B. Role of interlayer interactions 

To explain this, we must take into account the one 
remaining interaction present in the standard but not 
ideal model - the exchange J'^ between triangular planes. 
Like the (now absent) D term, it is strongly relevant, 
and unlike the D term, it is not averaged out in any 
field orientation. The appropriate continuum limit has 
already been given, the dominant piece being the last 
term of Eq. (57), reproduced here: 



(71) 

with -f1 = J" /2. Written explicitly using bosonization, 
it becomes 



(72) 



with 7^ = 2A\^'j^. This coupling is more relevant than 
either the SDW or cone interaction, in the entire range 
of magnetization. We may therefore expect that it scales 



to strong coupling unambiguously before any competing 
interactions. To check this, we again consult the com- 
parison of critical temperatures in Fig. [Sj Neglecting the 
effects of the D term, as we have just discussed, we see 
that the interlayer interaction is clearly dominant for all 
magnetizations below about 80 percent of the saturation 
value. Above this magnetization, the larger bare value 
of the cone interaction, which is of order J' rather than 
J" ^ overcomes the difference in scaling dimensions and 
controls the physics. In this high magnetization regime, 
the physics is therefore very similar to that described in 
the previous sections, and an incommensurate cone state 
is expected. 

In the remainder of this section, we focus on the main 
field regime, where 7^^ is dominant. The latter obeys 
the RG equation (using the dimensionless coupling 7^^ = 
^'j^/{vK\) as discussed in Sec. II D) 



da'i = (2 - 2A±)7: 



(73) 



where A± = iiR^ is the scaling dimension of the 
fields. Integrating this to the scale ^ = ln(^'') such that 
^±{() ~ V defines the length scale 



(74) 



For lengths shorter than one-dimensional fiuctuations 
are significant and approximately those of free chains. On 
longer length scales, we expect that 7^^ drives ordering 
of the 5^2;7r fields. H2 in Eq. (71) is minimized by con- 
figurations of the form 



(75) 



where |?/^| is a real number giving the magnitude of the 
spontaneous moment, and ii}y is a classical phase that 
can be chosen independently for each vertical a-b plane 
specified by y. Note that longitudinal order is strongly 
suppressed at this scale, (Sy^^) = 0, by the uncertainty 
principle (in bosonization this follows from the duality 
of the and (p fields). We expect by scaling that ~ 
^^//^-A± ^ -^^ refiecting the suppressed magnitude of 
magnetic order by Id fiuctuations. Hence 



IVI = <t(M) 



(v) 



2-2A_|_ 



(76) 



where t he p refactor (j{M) is computed by CMFT in Ap- 
pendix D5 We estimate ^ 0.25 — 0.3 over most of 
the field range. From this, we can estimate the lowering 
of the energy density due to the establishment of such 
two-dimensional order, simply by taking the expectation 
value of Eq. fnl : 



(77) 



An experimental measure of this energy density is the 
XY spin stiffness along the a axis, which is of the same 
order, pg-a ^ ^2d- Note that the spin stiffness along the 
b axis is much larger, of order pg-^b ^ 
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C. Exchange coupling of a-b planes 

The arbitrary choice of i^y for every y is a consequence 
of the fact that the dominant interaction, 7^^, does not 
couple different a-b planes. Less relevant interactions can 
and do remove this arbitrariness, ultimately determining 
the precise nature of the ordered state. 

To study this, we first include exchange interactions 
between chains within the b-c planes. In the standard 
model, this is only the J' coupling along the nearest- 
neighbor diagonals. However, it was argued in Ref. |14| 
that it is important to also take into account weak ex- 
change J2 between spins on second-neighbor chains sep- 
arated by distance Ay = 2. While clearly J2 <C J^ it is 
important because it is unfrustrated, unlike the J' inter- 
action. 

In the continuum limit, these couplings lead to the 
Hamiltonian 



y,z 

+725+^;^5-+2,^;7r+h.C.}, 



h.C. 



(78) 



with 72 = J2/2 (and 7cone = -^V^ as given earlier). 
Taking the expectation values using Eq. (75), the re- 



sulting renormalized Hamiltonian can then be treated 
classically, and minimized to find the ground state. It 
is evident that the "twist" interaction 7cone favors an 
incommensurate state with kx ^ tt. To describe this re- 
quires allowing for non-zero gradients dx'^y While such 
configurations are not ground states in the absence of 
7cone5 they are low in energy, because a small gradient 
comprises a soft (Goldstone) mode. The magnitude of 
the associated incommensurability is determined by a 
balance of 7cone with the gradient terms in i^o, which 
of course favor commensurate order at /c^^ = tt. We 
therefore include the latter, and write the entire effective 
Hamiltonian explicitly in the i^y variables, which we will 
allow to be x-dependent but independent of z according 
to Eq. ([75|). The total energy becomes 



y 

-gcone{dx^y ^ d^l^y^i) COs[/3{l9y - ^^^^+1 )] | , (79) 
with g2 = 272|V^P and ^cone = TconelV^P/^- 

Now we can see that, for ^2 > 0, which is expected 
from antiferromagnetic superexchange, the two interac- 
tions strongly compete. In this case, the minima of the 
g2 term are states with 



Try 6mod(jy,2) 

2^ ^ ' 



(80) 



Here 9o, 61 define the overall phase on the even and odd 
chains, respectively. Inserting this into the twist term, 
one finds a vanishing result due to cancellations when 



the sum over y is carried out, even if 69 and 61 are 
allowed to have gradients. Hence this solution has en- 
ergy density equal to —g2- This is a commensurate "an- 
tiferromagnetic" (AF) state. Conversely, the solutions 
which minimize the twist term have "dy = tzx, for which 
the g2 term is maximized rather than minimized. Here 
= 2gcone/^ IS determined by minimizing the full energy, 
leading to the energy density — 2^cone/'^ + ^2- This is the 
incommensurate cone state. Comparing the energies of 
the two states, one finds that the AF state obtains for 
92 > 9 cone This requires a minimum value of second 
neighbor exchange for the commensurate state, J2 > J|, 
where 



T* _ 

J2 — 



(81) 



For CS2CUCI4, J| is very small, and is in fact only a 
few percent (< 5%) of J for the relevant field range. 
Moreover, we argue below that the above value of J| 
is actually an overestimate, as it neglects a fluctuation- 
generated interaction which is of the same order. Thus, 
an exceedingly tiny second neighbor coupling J2, likely 
undetectable directly, qualitatively changes the ground 
state. In general, re-expressing the minimum energy den- 
sity in terms of bare variables, we have 



£3d--Max|j2|V^P, ^ 



n 



(82) 



This energy scale determines the spin stiffness along the 
c axis, ps-c ^ ssd- 

D. Locking of even and odd a-b layers 

When J2 is dominant in establishing three-dimensional 
AF order, it is ineffective in coupling the even and odd 
layers. As a consequence, there remains an artificial de- 
generacy of solutions, specifically, one may make opposite 
rotations of the phases Oi and 02- This rotation is not 
a true symmetry of the microscopic theory. However, 
the simplest possible coupling of phases in neighboring 
chains, of the form cos P {i^y — ^y-\-i)^ is prohibited by 
reflection symmetry, see (B14) and (B17). Instead, the 



leading possible coupling between neighboring chains is 
of the form 



^bq = +^bq Y jdx COs[2/3(^^,^ 



^.+1,.)]- (83) 



Here we have already assumed Eq. (75), and taken the 



average of the fluctuation-generated interaction. See Ap- 
pendix |E] for details. For classical XY spins with phase 
this interaction would correspond to a biquadratic 
coupling {Si • Sj)'^^ between spins on neighboring chains. 
Such fluctuation-generated biquadratic interactions are 
indeed familiar from the theory of frustrated magnets, 
and are a manifestation of "order by disorder" .^^ In 
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that context, it is well-known that fluctuations generally 
favor collinear states, which requires ^bq > 0. This is 
indeed conflrmed by the microscopic calculation in Ap- 
pendix IE) which leads to the estimate 



^bq 



i\2 



(84) 



In the AF phase, we may use the solutions for Oy^z 
determined above, and hence rewrite Eq. ([83| as 

(85) 



EhG, = -ghqLxLyLz cos[2(6o - 6i) 



Clearly states with 0o = 0i + nir are preferred, which 
implies commensurate, collinear, AF order. The conden- 
sation energy density associated with the selection of the 
collinear order is thus ^cit ~ ^bq (the reason for the 
choice of this subscript will become clear in the next sec- 
tion). Physically, this energy scale determines the gap of 
the antisymmetric pseudo-Goldstone mode correspond- 
ing to Oo — Oi, which is of order 



(86) 



This can potentially be measured as an emergent low 
energy (but gapped) mode in neutron scattering. 



Comparing Eq. (82) with Eq. (84) and Eq. (85), we 



observe that the energy gain due to ^bq term is of the 
same order as the energy gain of the incommensurate 
cone state. This is not a coincidence as the both effects 
have their common origin in the inter-chain exchange J^ 
This suggests that, even in the absence of any micrscopic 
J2 exchange, a collinear state could be energetically pre- 
ferred to the cone state. However, the RG approach used 
to obtain Eq. (84) is not accurate in determining the 0(1) 
numerical prefactor, which is essential for making such a 
comparison quantitatively. Thus we at present can only 
speculate that this might be the case. Even if not, these 
considerations imply that the interaction J2 needed to 
induce the AF state is even lower than the estimate in 



Eq. (81). 



VI. FIELD ALONG c AXIS 

Experimentally, this fleld orientation shows the most 
complex phase diagram. In addition to the commensu- 
rate "AF" phase seen for flelds along the b axis, a broad 
region of incommensurate phase is also clearly observed 
in NMR measurements^^ (and defined by earlier mag- 
netization measurements^). Within the model used up 
to now, the difference in phase diagrams for fields along 
the b and c axes is inexplicable: the Hamiltonian has a 
symmetry under spin rotations within the b-c plane. 



interactions in Sec. Ill Cl for our consideration. As we saw 
in the previous section, DM terms whose D-vector is or- 
thogonal to the applied field average out rapidly in the 
presence of an applied field. Hence we need consider only 
components of the D- vectors along the c axis. There are 
two independent such terms: Dc and D^^. Given that 
DM terms are generally proportional to the correspond- 
ing exchange, we expect Dc to be the largest of the two, 
and we focus on its effects (in Appendix F3 we explain 
in detail why D^^ can be neglected). It introduces the 
perturbation 

He - I)c^(-l)^ jdxJl,{x) 

y,z 

^ 4^(-l)^ jdxdjy^z, (87) 

where dc = vDc/{l3J). 

Notably, He is linear in the boson fields, and hence, in 
the absence of any other interactions, the term dc can be 
taken into account exactly. Moreover, it is actually a pure 
boundary term, whose effect on the energy depends solely 
on the winding numbers^ [Oy^z{^) — ^y,z{—^)]/{'^^P)^ 
and vanishes in the zero winding number sector. How- 
ever, He favors sectors with non- vanishing winding num- 
bers (proportional to Lx, in fact). 



B. Dc does not compete with J,J\ and J'' 

To understand the degree of competition of He with 
the other interactions, it is instructive to consider the 
shifted variables 

eyAx) = 9y,,{x) + {-l)y^x. (88) 

With this shift, dc is "eliminated" from the free Hamil- 
tonian, up to a constant: Ho[0] -\- Hc[0] = Ho[0] + const.. 
Physically, this change of variables corresponds to a 
shift of the dominant wavevector of correlations from 
kx = 7T to kx = TT ^ Pdc/v = TT ± Dc/ J- Significantly, 
the dominant 7^ coupling is invariant under the shift: 
H!2[0] = H^[e]. Thus J'' and Dc do not compete. The 
same is true for the 72 (J2) interaction. Thus, in the 
region where the AF phase appears for the case of field 
along the b axis, since the J^' and J2 couplings domi- 
nate, we expect the energetics is unchanged at the three 
highest energy scales in Fig. [9] 



C. Commensurate-Incommensurate Transition 



A. DM interaction on chain bonds 

Therefore additional spin-rotational symmetry break- 
ing interactions must be included to explain this discrep- 
ancy. We therefore turn to the general set of allowed DM 



Differences do appear, however, once the cone and bi- 
quadratic interactions are considered, as these are not 
invariant under the shift in Eq. (88). We focus on the 



putative AF region, for which we may assume the de- 
composition in Eq. (80). Allowing for small gradients 
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in 9±(x, 7/, z) = 60 (a:, 2;) =b ©i(x, 7/, z), we obtain the 
continuum hamiltonian H = Hj^ + H- , with 



(89) 



-ffbq COS 



(2e_)}, 



(90) 



where Cx = '^/4/3^, Cy = and = 7^^/4. 

Here the cone interaction has dropped out, and the 
low energy Hamiltonian has decomposed into two de- 
coupled parts. The first, is simply the Hamilto- 
nian of a free massless boson. It describes the Gold- 
stone mode 6+ associated with spin rotations about the 
field axis. The second part, is the familiar sine- 
Gordon model, discussed earlier in Sec. |HID4l In this 



case it is in three dimensions, but this has no signifi- 
cant consequences. As in Sec. HID 4[ the sine- Gordon 
model describes a commensurate phase (here, the AF 
state) and an incommensurate (IC) one, separated by a 
Commensurate-Incommensurate Transition, or CIT. 

The results for the CIT can be taken over directly from 
Sec. |III D with the mapping 0_ (p/2. Here, it is dc 



which plays the role of the tuning parameter, favoring 
the IC phase, for \dc\ > <i*, where 



On entering the IC phase, the system forms a soliton lat- 
tice, with a corresponding incommensurate wavevector 
qo (measured relative to the AF state). Note that the 
IC phase found here is thus a smooth deformation of the 
AF state, which makes it quite distinct from the cone 
state, which is also incommensurate. The incommensu- 
rate wavevector grows rapidly after the CIT, which can 
be seen by translating Eq. (56) to the current case: 



^0 



1 



ln[(|4|-4*)R*] 



(92) 



where the brackets (dxO-) indicates the spatial average. 
Once the above logarithm is not large, the solitons are 
strongly overlapping, and Eq. (92) is no longer valid. In- 



stead, one may simply minimize the energy neglecting 
^bq, which gives 



Pd, 

V 



Dc 
J 



(93) 



To summarize, varies from its maximal value given 



in Eq. (93) at the low field end of the IC phase, and 



decreases with increasing field, vanishing asymptotically 
according to Eq. (92) at the CIT to the AF phase. Be- 



cause the variation in Eq. (92) is so rapid, very likely qo 



appears approximately constant in most of the IC phase, 
dropping precipitously to zero in a narrow region near 
the CIT. 



VII. EXPERIMENTAL CONSEQUENCES 

In this section, we consider a few key experimental con- 
sequences of the analysis of the previous sections. First, 
we give explicit expressions for the spin structures in the 
various phases predicted there, which should be useful for 
comparison to neutron scattering measurements. Next, 
we derive the nuclear magnetic resonance (NMR) line- 
shapes in each of these phases, using these expressions, 
and compare to experiments by Takigawa and collabo- 
rators. Finally, we describe the phase diagrams in the 
magnetic field-temperature plane, for the different field 
orientations. 



A. Explicit spin structures 

Here we reconstruct explicit formulae and plots of the 
spin ordering patterns in the various phases discussed 
earlier. 



1. Cone state 

First consider the incommensurate ordered "cone" 
state, described in Sec. |IVB[ which occurs for any field 
along the a axis. This is described by Eq. (66). Using 



Eq. (24), we can express the microscopic spin operator's 



(91) expectation value 



(94) 



= {-i)y'N cos[(^ + {-iyqo)x + e,] 

= {-l)y'Ncos[i7r + qo)x + ((-1)^ - l)7rx + (-1)^6,]. 

To further simplify, we note that the x coordinate takes 
integer values for even y and half-integer values for odd 
As a consequence, ((— 1)^ — l)7rx is an integer multiple 
of 27r unless y is odd and z is odd. This allows this factor 
inside the cosine to be removed in favor of an overall 
{ — factor in front of it, which cancels the one already 
present. Therefore one finds, finally 



\ -i^yz / cone 



7Vcos[(7r + qo)x + 9;^], 



(95) 



where O;^ = { — l)^Sz' Similar manipulations for the y 
component of the spins give 

{Sly.),,,, = {-irNsm[{7r + qo)x + e,]. (96) 

and of course, one has 



I cone 



M. 



(97) 



2. Antiferromagnetic phase 



Here we simply apply Eq. ( 80 ) , and use Oq 
O + 7rn (with n = 0, 1) as preferred in the commensurate 



0,01 
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"antiferromagnetic" (AF) phase by Eq. ( [85| ). This gives 

\^xyz/ X¥ 



(-l)^IV^I cos[7rx + cT7Ty/2 + 6], (98) 
(-l)^IV^I sin[7rx + a7ry/2 + 6], (99) 



where the a = ±1 = (—1)"^. These equations describe 
a state in which the x — y components of the spins are 
coUinear. This may not be obvious, but is true because 
the combination 7rx-\-a7ry /2 is always an integer multiple 
of TT, owing to the fact that x is integer (half-integer) for 
even (odd) y. Combined with the constant uniform mag- 
netization, Eq. (97), these equations describe a co-planar 



spin state, distinct from the three-dimensional cone con- 
figuration. We note, however, that small perturbations 
due to the various DM interactions will probably disrupt 
this ideal coplanarity. The commensurate nature of the 
ordering is, however, robust. 



simply proportional to the magnitude of the total effec- 
tive field. The simplest approximation, which we take 



here, is to assume in addition that \hf \ <C |h^^*|. This is 
certainly so in intermediate and high-field regions, which 
we focus on. For lower fields, of the order of 1 — 2 T, this 
may not be such a good approximation.^^ But even in 



this case the off-diagonal contribution (see (105) below). 



which is central to our consideration, should be smaller 
than the diagonal one and a modified expansion in off- 
diagonal components of h-^^ should be possible. With 
these assumptions in mind, and disregarding the g-factor 
anisotropy (so that h^^* cx H), we can approximate the 
shift due to the hyperfine interaction by 



H, 



(103) 



where H = H/|H|. In CS2CUCI4, there are two inequiva- 



Incommensurate phase for fields along c axis 



Here we consider the incommensurate phase which is 



discussed in Sec. |Vl) For simplicity, we will ignore the 
narrow but subtle region in the vicinity of the CIT, where 
a non-trivial soliton lattice should be taken into account. 
The basic symmetry of this phase is well described by the 
"smooth" regime (corresponding to strongly overlapping 
solitons), where we simply treat the incommensuration 
as linear shift of the phase fields, i.e. we take Oy^z in 
Eq. (88) as constant. The preceding formula now are 
modified to 




FIG. 10: Scheme of the trasferred hyperfine interaction for 
Cs(A). The signs +/— re fer t o the relative signs of the off- 



diagonal entry nij in Eq. (105). nij takes equal value for two 



where = (3dc/v - see Eq. (93) - and Oy is a phase 



taking two distinct arbitrary values for even and odd y. 



sites of a triangle as shown by the dashed lines while that for 
y^j, (100)^^^ other site can be different as shown by the dotted line. 

{Sly,\c = (-l)'IV'l sinKTT + (-l)yqo)x + ny/2 + 9,], (101) 

lent Cs sites measured in the Cs NMR measurements. We 
focus on the Cs(A) site, which is located slightly above or 
below the center of a triangle of Cu spins. The NMR re- 
sponse of the^^ Cs(B) site, which is coupled to Cu spins in 
the two planes adjacent to it, is more difficult to analyze 
at present, but we do expect the qualitative conclusions 
derived below to remain valid for this situation as well. 
In the Cs(A) case, the hyperfine field of nucleus i located 
between spins j at the sites of this triangle should be 
given by the sum of three transferred contributions: 



B. NMR lineshape 

Recent NMR experiments by Takigawa and 
collaborators^^ have revealed numerous phases and 
transitions in CS2CUCI4 in magnetic fields. Here we wish 
to address the signatures of the phases predicted in this 
paper in the NMR lineshape. The basic approach is to 
consider the Hamiltonian of a given nuclear spin to 
be the sum of two effective fields 



(104) 



J nn z 



Ha = (hf* + h^^) • I 



(102) 



where h^^* is the effective field on the nucleus i due to 
the external field H, factoring in any anisotropics of the 
nuclear g-tensor (which are believed to be smal^). The 
remaining "hyperfine field" h^^ represents transferred hy- 
perfine interactions with nearby electronic spins. The 
NMR resonance frequency of this particular nucleus is 



where K^j is a tensor describing the anisotropic trans- 
ferred hyperfine exchange from the Cu spin at site j 
to the nucleus i (see Fig. 10). According to recent 



measurements,^^ in these tensors, the only significant off- 
diagonal entry is [K^j]^^ = [K^j]^^ = hZij^ and, moreover 
f<iij takes equal values for the two sites j of the triangle 
which are on the same chain. 

For magnetic fields along a and c, this off-diagonal 
transferred exchange is crucial in determining the NMR 
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lineshape. Let us see how this occurs. In either of these 
cases, we define, as usual the z axis of spin along the field 
axis. Let us then take the x axis of spin along the other of 
the two, i.e. for H = a, take = . and conversely, if 
H = c, then 6'^ = 6'^. From Eq. (Il03|, the NMR shift is 



Cone state 



entirely determined by the z component of the hyperfine 
field. This, in turn, is given by 



E 

j nn i 



(105) 



In all of the phases predicted for CS2CUCI4, the compo- 
nent of the spins parallel to the field is constant, and 
equal to the average magnetization M (this is not true in 
the SDW phase, which is expected in the ide al 2d case of 
Section [nil ). Therefore the first term in Eq. ( |1Q5[ ) gives a 
constant contribution to the shift, which is the same for 
all Cs(A) nuclei. Thus 



Aui (X const. + ^ K:ij{Sj). 



(106) 



Using the experimentally determined form of the hyper- 
fine couplings, and dropping the constant, one has 

(107) 

Here we have absorbed the proportionality constant in 
the shift into the definitions of a^i and i<i2- We are now in 
a position to evaluate the NMR lineshape for the different 
magnetic phases. 

a) Cone state 




Sv 



-A A 

b) AF state 



c) IC state 



Sv 




Sv 



FIG. 11: Schematic NMR spectra in (a) the cone state, (b) 
the AF state, and (c) the IC state. On approaching the com- 
mensurate AF phase from the IC phase, i.e., ^0 ^ limit, 
four peaks merge pairwise as indicated by the arrow. 



In the cone state, we can use Eq. (95) to evaluate 
Eq. (|1Q7[). One obtains 



{-i)yN [^2 - 



2f^i sin f j 



cos[(7r -^qo){x - 



(108) 
1/2) + 6,]. 



Now the NMR lineshape refiects the distribution of shifts, 
p(Az/), over all the Cs(A) sites. We may consider this as a 
sum of distributions of the shifts for the nuclei associated 
with each pair of chains, i.e. ranging over x for fixed y and 
z. Because qo is incommensurate, the argument of the 
cosine above is distributed uniformly over the full angular 
interval from to 27r. Thus the cosine itself is distributed 
between —1 and +1, and we obtain a distribution for the 
shift, for fixed y and z with support between ±N\hZ2 — 
2ki sinf |: 



-,e[A-\Au[ 



with 



N\k2 



2ki sin^l. 



(109) 



(110) 



We see that the distribution is in fact independent of y 
and z^ so that the full distribution over all Cs(A) sites 
is identical to that for a single pair of chains. It has 
two peaks, at the edges of the distribution, = ±A as 
shown in Fig. 11 ^a). 



2. AF state 



Applying Eq. (p8| to Eq. (107), we obtain 



Au, 



cr(-l)^+^|V^|A^2 cos[7r(x + (jy/2) + 6] 

±A^|V^|cOSe. (Ill) 



One expects therefore two sharp peaks in the Cs(A) NMR 
spectrum, separated by 2/^2'0|cos6| (see Fig. [iTJb)). 
Note that cos 9 is generically non-zero, as argued by sym- 
metry in Appendix |C 2[ 



3. IC state 



Here we apply Eq. (|100l) to Eq. (|107|). We find that 

(112) 



the NMR shift can be written 
where 



2 qo 



IV^I 4^1 sin 
-{-l)y 41^11^2 sin f sin(6o + 61) 



(113) 
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and 



tan = 



(114) 



/^i(sin(go+Bo)-sin Oo)+/^2 sin(^-ei) 
/^i(cos(qo+0o) — cos Oo)+/^2 cos(%— Oi) 
/ti (sin(go — 0i)+sin 0i)+/^2 sin(^+Oo) , , 

ub y odd 

Ki (cos(go — ©i) — COS Oi)+h:2 cos(^+Oo) 



y even 



For each we expect from Eq. (112) a continuum hne- 



shape of the form of Eq. (109), owing to the incommen- 
surate wavevector qq. However, in general. Ay takes two 
distinct value s for even and odd y (owing to the ( — 1)^ 
factor in Eq. ( |113| )). Note that the prefactor of this term 
is non- vanishing since si n(8Q + 0i) is generally non-zero, 
as argued in Appendix |C 3[ Hence the full Cs(A) line- 
shape is expected to be the sum of both distributions, and 
hence has four peaics, at Au = ±Ao,±Ai as schemati- 
cally shown in Fig. 11 'c). Notably, these peaks merge 
pairwise as qo 0, i.e. on approaching the commensu- 
rate AF phase from the IC state. Precisely such a merg- 
ing of the peaks has been seen in the NMR experiments 
by Takigawa and collaborators^ 



However, we note that the nature of the phase transi- 
tion is somewhat subtle, and probably not properly de- 
scribed by this approximation. Neglecting J" ^ which has 
a very weak effect upon the cone state (see Eq. (69)), 
the system is effectively two-dimensional, and as a con- 
sequence exhibits strong effects of thermal fluctuations. 
Since the DM terms (e.g. Dc) with DM-vectors perpen- 
dicular to a are also negligible here, the Hamiltonan has 
approximate XY spin rotation symmetry. As a conse- 
quence, the cone phase is approximately a quasi-long- 
range-ordered state at T > 0, and its thermal transition 
should be of Kosterlitz-Thouless (KT) type. Obviously 
the CMFT approximation does not describe the KT uni- 
versality class, and instead predicts mean-field critical 
behavior. 

It is interesting to verify that, nevertheless, the magni- 
tude of Tc obtained from CMFT agrees with an analysis 
based on KT theory. For simplicity, we will focus on the 
DM-dominated field range, and neglect entirely J' and 
J^' for simplicity. In this case, the system decouples into 
2d triangular x — y planes, consisting of chains connected 
by the DM interaction only. Taking the expectation value 
using Eq. (59) in Eq. (58), one obtains the energy 



C. Phase diagrams 



In Sees. IV|V|VI[ we have determined (most of) the 



zero temperature phases for the three major field orien- 
tations. Here we discuss the extension of these results to 
T > 0. 



H2d = Y.I ^^{^(S:.^y,.)'-21^l^^l'cos/3(^^,,-^^+i,,)}. 

(115) 

In the ordered phase, one may expand the cosine, and 
take the continuum limit (in y) for fields i^y^z that are 
slowly- varying in y: 



1. Field along a axis 

This is the simplest case. At zero temperature, the 
cone state extends across the entire field range from zero 
up to saturation. We have seen that it is predominantly 
controlled by the DM interaction D = D'^^ perturbed 
somewhat by the interchain exchange J' . We therefore 
expect a single phase boundary, Tcone{H). One estimate 
for this curve is obtained from CMFT, and is shown 



in Fig. 19 (plotted versus magnetization M rather than 
field). One observes that Tcone at first increases with the 
applied field for small fields, and then reaches a maxi- 
mum, followed by a decrease to zero at the saturation 
field. 

These trends can be understood simply as follows. 
With increasing magnetization, the spins become more 
XY-like, which decreases the scaling dimension A±. As a 
consequence, the DM interaction becomes more relevant 
with increasing field, enhancing the critical temperature. 
However, on approaching saturation, the magnitude of 
the transverse components of the spins, which constitute 
the cone order, decrease to zero, and hence suppress the 
ordering temperature to zero. 

As these trends are correctly captured by CMFT, we 
may perhaps trust the result for the phase boundary. 



2d 



I dxdy {"^{d^^^f + D\N\^pHdy^,f}. 



(116) 

Now, according to KT theory, the critical temperature 
is proportional to the geometric mean of the two stiff- 
nesses, i.e. Tkt ^ VvD\N\. We are neglecting all 0(1) 
prefactor s here, as we are only interested in the scaling 
behavior. Now fro m sc aling, or from the CMFT calcula- 
tions in Appendix 



D5 



one has \N\ - {D/v)^^/^^-^^^\ 

from which one obtains finally Tkt ^ v{D /vY^^'^~'^^^\ 
Precisely the same scaling is found directly from the 
CMFT treatment at T > in the Appendix [dI] (Please 
note that the described calculation corresponds to setting 
J' = the coupling Fcone^ ( D27 ). That implies (70 = 
which, via second equation in ( |D23 ), leads to the scal- 
ing Tc ^ v{D /vY^^'^~'^^^\) It may appear surprising 
that the two approaches, which describe the transition 
so differently, agree in this respect. The reason for the 
agreement is that the scale of Tc is entirely determined 
by the scaling properties of the weakly perturbed one- 
dimensional chains. Any approximation which respects 
this scaling (and both the CMFT and the KT analysis 
do) will obtain the same order of magnitude answer. Dif- 
ferences would appear in the prefactor, which, however, 
is beyond the scope of the rough KT analysis carried out 
here. 
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2. Field along b axis 

The situation in a field along b is considerably more 
complicated. The analysis in Sec. [V| implies at least four 
phases at T = 0: the zero field phase "spiral" phase, dom- 
inated by the AF phase, the high-field cone phase, and 
the saturated phase. Due to the difficulty of treating the 
competition between the D term and magnetic field when 
the two are comparable, the intervening range between 
the AF and zero field phases has not been fully clarified 
here. Thermodynamic measurements^ appear to show a 
single transition between the spiral and AF states, and 
therefore the absence of any intermediate states. 

For the "high field" phases (i.e. in the region where the 
D term is negligible) , we can attempt to apply CMFT to 
determine the uppermost phase boundaries, describing 
the transitions from the ordered to paramagnetic states. 
The key observation is that both ordered phases ( AF and 
cone) are driven by the same J'^ interaction. The two 
states are only distinguished by the competing effects of 
the weaker (at least in the renormalized sense) J' and 
J2 interactions. Thus the upper phase boundary should 
be approximately continuous across this field range up to 
saturation, and not very sensitive to the precise nature 
(AF or cone) of the ordered phase it demarcates. This 
boundary should be similar in shape to the Tcone{H) dis- 
cussed above, as it arises from a term of the same scaling 
dimension as D in that case, and suffers the same reduc- 
tion on approaching saturation. 

The high-field region requires one further phase bound- 
ary, between the AF and cone states. This should be 
approximately vertical, but is expected to bend "to the 
left", as the cone state has higher entropy than the AF 
one, and is thus favored with increasing temperature. 
This reasoning is based on higher order effects of other 
DM interactions that we have negl ected up to now (men- 
tioned in passing in Sec. VII A 2[ ): these sub-dominant 
terms are expected to break the U{1) rotational symme- 
try of the AF state and gap out its Goldstone modes, 
while preserving the commensurability of the AF struc- 
ture. The cone state, being incommensurate, is expected 
to not be affected by these small perturbations and pre- 
serve its gapless excitations. As a result, we expect the 
entropy of the cone state to be greater than that of the 
AF one, and result in the mentioned bending of the AF- 
cone boundary to the left. The transition between the 
two states is first order, and observables such as the or- 
dering wavevector jump at the critical field. 



3. Field along c axis 

In this, most complex field orientation, all the phases 
predicted for the field along b must appear, and in addi- 
tion the IC state, taking up some territory betwee n the 
AF and spiral phases. Experiments seem to 
even beyond these 5 states, one or two additional ones 
in the regime when the D term is comparable to the 



Zeeman energy. At present we have little to say about 
these states. The observed linear relation between the or- 
dering momentum of the intermediate "S" state (in the 
notation of Ref.'S) is suggestive of an SDW phase, but at 
the present we do not have a good understanding of how 
the competition between the DM and the Zeeman terms 
may bring out the SDW order discussed in Sec. |IIIB[ 

In the high-field region, we expect, by similar argu- 
ments to the previous section, a rather continuous bound- 
ary between the IC, AF, and cone states and the param- 
agnetic phase. The AF-cone boundary should appear 
very similar to that for this field orientation as well. 
The new feature introduced here is the IC-AF boundary, 
which is the location of the CIT. Like the cone phase, 
the IC phase is expected to have more entropy than the 
AF state, and hence be stabilized by increasing tempera- 
ture. Experimentally, this boundary bends quite sharply 
"to the right", in marked contrast to the nearly vertical 
AF-cone transition line. 

To understand this, consider the expression for the CIT 
location, Eq. (91). It can be rewritten, using the ex- 



pressions for ^ bq? Eq. ( 84 ) , and for dc (given in the text 
following Eq. (87)), as 



I3JJ'' 



(117) 



where we drop all 0(1) factors. A simple treatment, 
which takes into account some of the one-dimensional 
fiuctuations, but not the high-dimensional ones, is to sim- 
ply apply Eq. (91), but replacing with its reduced 



value atT = Tcit- In general, this is difficult to actually 
calculate analytically, e.g. with CMFT, but the detailed 
form is not important to our argument. Scaling implies 
that it can be written 



(118) 



where Tc is the critical temperature at which ip vanishes, 
and the scaling function T is smooth and obeys J^ifi) = 1 
and J^{1) = 0. For illustrative purposes, we can take the 



simple approximation T{t) = 1 
one finds 



t . Using this form. 



^CIT 



vDr 



(119) 



By construction, the right-hand side vanishes at the zero 
temperature CIT, where Tcit = 0. Both ipo and /3 are 
rather weak functions of magnetic field. However, the ve- 
locity V varies considerably (on the scale of J) with field, 
indeed vanishing as saturation is approached. Hence, 
the right-hand-side increases rather quickly with field, 
leading to rapid variation of Tcit with an approximate 
square-root form, consistent with experiments. 

It is interesting to note that, in experiment, §^ the 
AF-IC and AF-cone boundaries are observed to approach 
each other very closely with increasing temperature, lead- 
ing to an extremely narrow range of transition directly 
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from the paramagnet to the AF phase. This suggests 
some physical mechanism which "avoids" this transition. 
In fact, one can argue that, according to Landau theory, a 
continuous AF-paramagnetic transition is forbidden for 
this field orientation. To do so, consider the Landau ex- 
pansion of the free energy F in the "order parameters" 
ipy^zi^) = {'^y,zi^))- Such an expansion, in powers of 
ipy,z^ is valid near any putative continuous transition. We 
presume i/jy^zi^) to be a slowly- varying function of x. The 
Landau expansion has the form F = F2 + F4 + • • • , where 



c.c. 



h.c. 



-idc{-iy^ljl^zdx^y,z^l2 {i^lz%+2,z 

+7± Kz^y^z^i + h.c.) + r\^y,zf} (120) 
contains quadratic terms in the order parameter, and 



E 

y,z 



J dx^ultpy^. 



+7bq[(V^;,.Vl,.)'+H } 



is quartic. In Eqs. ( 12Q|121 ), the couplings with overlines 
on them are analogous to the corresponding couplings in 
the bosonized Hamiltonian, as can be seen if one assumes 
i^y,z = iV^le*^^'^. They are, however, from the present 
point of view, phenomenological coefficients which are at 
best proportional to those microscopic couplings. 

Let us consider possible continuous transitions from 
the paramagnetic state. In this case, we may assume 
\i^y,z\ is arbitrarily small, and thus F4 is a small per- 
turbation to F2. The transition occurs on decreasing 
r from large positive values, at the point at which the 
smallest eigenvalue of the quadratic form in F2 van- 
ishes. Fourier transforming into the two-site unit cell, 
i^yA^) = Jd^k/{27rfilja,ke'^-''^'^yy^'^^^ with a = for 
y even and a = 1 for y odd, one obtains 



d^k 
(2^ 



K,k^ab{k)tlJb,k, 



(122) 



where the matrix T{k) can be decomposed into the iden- 
tity matrix, J, and the Pauli matrices, cr^, according to 



with 



^0 = 2 + COS kz + 272 COS 2ky + r. 



(123) 



J~ X — 2'y ckx COS ky^ 
T — d h 



(124) 



One immediately concludes that the minimum eigenvalue 
of T is 



^min = ^0 — ^ ^z 



This, in turn, should be minimized over k. Minimization 
over kx and kz is simple: the minimum occurs at 



\kJ = 



and kz = 7T. Then 



(126) 



^ . , N ,1 cPr 27^ cos^ ky 

J^miniky) = r-2il-^-^^ ^ + 272Cos2A:j, 

cP + 2^2 ^2 
= r-2ii- ^^^^ + (272 - y ) cos 2ky. 

(127) 

From here we immediately see that the minimum free 
energy is obtained for ky = when ^^jv > 272 and for 



7r/2 otherwise. In either case, we see from Eq. (126) 



that kx 0. The two cases therefore correspond to the 
cone and IC states, respectively. Noteably, the commen- 
surate AF state is absent. This is easily understood since 
it is stabilized by the biquadratic coupling ^bq, which 
in Landau theory corresponds to the quartic interaction 
7bq. Since this becomes parametrically small relative to 
the quadratic terms as |?/^| ^ 0, it cannot stabilize a com- 
mensurate phase in this limit. Thus we conclude that if 
there is a continuous transition from the paramagnet to 
an ordered state for this field orientation, it can only be 
to the cone or IC phases, and not to the AF state. Con- 
versely, if there is a direct transition between the para- 
magnet and AF states, it must be first order. This lat- 
ter scenario appears to be the case in experiment.^ We 
note that for fields along the b axis, where dc = 0^ di direct 
continuous transition to the AF state is possib le, si nce in 
that case ky = 7r/2 and kx vanishes from Eq. (126). 



Mk) - I IV cos2 ky^dl (125) 



D. Cs2CuBr4 

It is instructive to compare the case of CS2CUCI4 ex- 
tensively reviewed here with that of its isostructural 
equivalent Cs2CuBr4. The latter material is more two- 
dimensional, with J' / J ^ 0.75 as estimated in Ref. 50 by 
comparing the observed momentum of magnetic Bragg 
refiections, qo = (0, 0.575, 0), with the result of the series 
expansion calculations in Ref. EU This estimate should 
be taken with some caution, as the theory neglects DM 
coupling, which is clearly present in experiment (as wit- 
nessed by the distinct differences between the behavior 
in a field along a and perpendicular to it). Ref. [50] ar- 
gued that the inter-plane coupling in Cs2CuBr4 is weaker 
than in CS2CUCI4, since the ratio of the saturation field to 
the Neel temperature is approximately 1.5 times larger 
in Cs2CuBr4 than in CS2CUCI4, and usually the Neel 
temperature in quasi-2d systems is expected to be de- 
termined by inter-plane coupling. In principle, this need 
not be the case when DM interactions are strong, but we 
believe the conclusion is probably correct. Thus, com- 
pared to CS2CUCI4, we surmise that J' / J is increased 
and J" / J is decreased. This behavior is in agreement 
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with the estimate based on the band- structure calcula- 
tion of material's microscopic parameters in Ref. |52] 

We believe that reduced three-dimensional coupling is 
the primary reason for the observed cascade of phase 
transiti ons i n Cs2CuBr4 subject to magnetic field in b- 
c planePSl^^H^ Particularly striking is the observation of 
a robust M = ^Mgat magnetization plateau as well as 
a hint of possible second plateau, at or near 2/3 of the 
saturation magnetization. 

In the quasi-one-dimensional approach adopted here, 
as discussed in Section IIIIDl the existence of the SDW 



state is a necessary condition for the plateau. Given 
that inter-plane exchange J^' strongly favors cone state 
over the SDW one, we understand that Cs2CuBr4 with 
its small inter-plane coupling is indeed a good candi- 
date for the magnetization plateau. One must remember 
that this argument is based on one-dimensional reason- 
ing, the validity of which in Cs2CuBr4 {J' / J = 0.75) is 
questionable. However, the final outcome of this - that 
the magnetization plateau is stable in the full range of 
J' / J <: 1 ra tio - is completely consistent with two recent 
studie^^^^ which approached the problem as a spatially- 
deformed two-dimensional one. 

It is worth noting that abrupt variation of the SDW 
ordering momentum Q on approaching the plateau value, 
described by Eq. (56), can be clearly seen in Figure 9 of 
Ref. |50l While this strong feature was interpreted there 
as an indication of a first order transition, our theory 
would predict very similar behavior from a continuous 
two-dimensional C-IC transition. 

It is interesting to contrast the physical scenario emerg- 
ing from the quasi-one-dimensional approach to what is 
expected based on semi-classical physics and the more 
isotropic regime.l^ Notably, the phases immediately bor- 
dering the 1/3 magnetization plateau in the latter case 
are not collinear SDW states. Instead, Ref. 57 finds com- 
mensurate planar or non-coplanar incommensurate dis- 
torted umbrella states. These states are connected to 
the uud plateau state by continuous phase transitions 
which however can be driven first order by residual DM 
interactions,'^ which generally allow for cubic terms in 
free energy expansion. NMR measurements^^ find that 
the states below and above the plateau are incommen- 
surate, but cannot distinguish SDW from distorted um- 
brella states. These experiments and other^^ also find 
some hysteresis at the plateau edges, which suggests first 
order transitions there. As we have discussed, one ex- 
pects second order transitions in the SDW case, so this 
probably suggests that SDW state does not occur in 
Cs2CuBr4. This is also supported by the neutron scatter- 
ing experiments,'^ which observe a smooth evolution of 
the scattering intensity from zero field up to the plateau 
edge. 

Despite the evident absence of SDW physics in 
Cs2CuBr4, it is still interesting to consider the predic- 
tions of our theory for the plateaux themselves. Apart 
from the persistence of the 1/3 plateau to small J' / J ^ 
the most striking outcome of our theory is probably that 



the second "strongest" candidate plateau is not at 2/3 of 
saturation but at 3/5 of it. This feature should be taken 
as another definite prediction of our work. 



VIII. DISCUSSION 
A. Resume 

In this paper, we have presented a fairly thorough anal- 
ysis of the low temperature phases of CS2CUCI4, obtained 
from a quasi-one-dimensional approach. The results ex- 
plain most of the specific heat, magnetization, NMR, and 
neutron data available. Several aspects are particularly 
remarkable. First, contrary to the popular view of this 
material as an "anisotropic triangular lattice antiferro- 
magnet", for magnetic fields within the XY plane, the 
strongest two-dimensional ground state correlations are 
within the a-h planes, perpendicular to the nominal tri- 
angular {b-c) layers! Second, in establishing the phase di- 
agram, we have argued for the critical importance of four 
different very weak interactions {D, J" ^ Dc and J2), only 
two of which have been generally recognized {D and J") 
in prior work. It is remarkable that interactions of a mag- 
nitude of only a few percent of the largest exchange can 
induce entirely new phases. Finally, we have discovered 
an heretofore unnoticed commensurate-incommensurate 
transition in this material, and located its telltale signa- 
ture in NMR experiments. 



B. Relation to previous work 

The subject of quantum antiferromagnetism on the 
triangular lattice is long and storied. Here, we will re- 
view various aspects of the problem discussed in the lit- 
erature which relate to this paper. Some of the earli- 
est work^^*^ applied the random phase approximation 
(RPA), using bosonization results for one-dimensional 
Heisenberg chains, to study the susceptibility in the para- 
magnetic phase, and estimate critical temperatures. In- 
deed, the RPA is equivalent to the CMFT used here, as 
far as predictions of the critical temperature are con- 
cerned, provided the same interactions are taken into 
account. At a more general level, this early work cor- 
rectly emphasized the importance of the one-dimensional 
regime. However, the analysis here (and in Ref. 14) is 
much more complete in a number of significant ways. It 
treats the ordered phases below Tc, takes full account 
of anisotropic DM couplings, and includes fiuctuation- 
generated interactions which are ignored within RPA. 
These effects rather dramatically alter the phase diagram 
of CS2CUCI4 from the expectations of Refs. 18 and 58l 

Much of the theoretical work motivated by CS2CUCI4 
focused on the inelastic neutron structure factor, address- 
ing experiments^ which observed very broad lineshapes 
and extracted dispersion relations for putative "magnon" 



peaks. Several groups applied spin-wave theor}^ 
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study the ground state (staggered) magnetization and 
the structure factor, including higher order corrections 
in 1/S. The low energy dispersion of the zero field 
magnon peak is reasonably well reproduced by this ap- 
proach, while higher energy portions are not. To fit 
them, requires "renormalizing" the exchange couplings 
by hand, in a manner inconsistent with other measure- 
ments (e.g. at high fields). Moreover, the large contin- 
uum scattering is not obtained in this approach. An- 
other series of works study the excitation spectrum of 
"magnons" using series expansions. I^^Jt^ In our opinion, 
because the ground state is ordered, and the series are 
constructed from such a starting point, they are fairly 
reliable in determining the energies of magnon-type ex- 
citations of the system. (Although they do miss impor- 
tant finite lifetime effects which can be quite large in 
non-collinear spin configurations.'^ Indeed, the results 
compare well to the dispersion of the peaks of intensity 
in experiment. This method provides a useful compu- 
tational tool, especially helpful in estimating exchange 
couplings. However, it does not elucidate the mecha- 
nism of magnetic ordering or provide a full description 
of all the excitations. Thus it is much less useful if the 
ground state is not known (as in much of the non-zero 
field experiments) , and it does not address the dominant 
continuum portion of the experimental spectra. Sev- 
eral theories approached the ex citation spectrum from 
more exotic perspectives ,'^2^^^^''^ postulating proxim- 
ity to quantum spin liquid phases and/or quantum crit- 
ical points. We believe there is little support for such 
proximate exotic phases from experiments on CS2CUCI4. 
Instead, the most compelling exp lanation for the neutron 
experiments comes from a iheoi^^^^^ in which the ex- 
cited states are constructed from superposition of a small 
number of elementary "spinon" excitations of the individ- 
ual Heisenberg chains. This approach quantitatively and 
qualitatively explains the main features of experiment, 
with no adjustable parameters. Its success is a strong 
argument in favor of the quasi-one-dimensional approach 
adopted here. 

Several works address the ground states of spatially 
anisotropic triangular antiferromagnets, and CS2CUCI4 
in particular. In zero magnetic field, the ideal prob- 
lem (d iscussed here in Sec. Ill) has been heavily 
studied.' ^ ' -^^ ' ^^'^^ ' ^-^ l '^-^ ' '^^' Many approaches find simply 
that, in the quasi- Id limit, the correlations between 
chains are extremely weak, and either regard this small 
J' / J region as a "spin liquid^' or are inconclusive as to 
the actual ground state.^^^^^^The most recent series ex- 
pansion calculations of Re f.l22l favor a spiral state, but do 
not make a definitive conclusion. The approach described 
here, applied to this case in Ref.[l4l predicts definitively a 
collinear ground state, arising from a rather subtle fourth 
order fiuctuation effect. Very recently, a numerical cou- 
pled cluster method^ also obtained this state. Unfor- 
tunately, because of the very weak fourth order nature 
of the stabilization of this phase, we expect this to be 
a somewhat academic result. The DM interactions in 



CS2CUCI4 (and likely in many other anisotropic triangu- 
lar systems) completely overwhelm the fiuctuation effect 



and result in a spiral state, as decribed in Sec. Ill 

In a non-ze ro applied field, there has been less 
effort fUME^^ Spin wave theory, applied to the "stan- 
dard" model of CS2CUCI4, has considerable success in 
reproducing many of the features of the ground state 
phase diagram. It does not, however, explain the 
broad regions of AF and IC phases appearing for fields in 
the h-c plane. An explanation of the region correspond- 
ing to the AF phase found here was given in Ref. l36l 
based on a dilute spin-fiip approximation. However, the 
state obtained there differs from our AF state. It is a 
non-coplanar commensurate state, with the spin compo- 
nents transverse to the field lying in orthogonal direc- 
tions on neighboring chains. Such a state would obtain 
for negative biquadratic coupling ^bq < in ([83|), and 
would probably have distinct signatures in NMR mea- 
surements. In our opinion, the sign ^bq > 0, found 
here, is rather more natural, and more consistent with 
the usual expectation that fluctuations ("order by dis- 
order") favor more collinear states.^-^ The IC phase and 
CIT discussed here are entirely new, and could not pos- 
sibly have been obtained in previous works, all of which 
assume rotational symmetry of the exchange interactions 
in the h-c plane. As a general point, it is not too sur- 
prising that spin- wave based approximations can capture 
most of the phases in an applied field. This is because 
we have found that these have "classical" order parame- 
ters (non-zero (<S^^(x))), and are thus adiabatically con- 
nected to mean-field states. Such approaches may, how- 
ever, have large quantitative errors, and furthermore, 
may miss states where fiuctuation-induced interactions 
are important (such as the AF state). The problems 
with spin-wave theory are most extreme actually in the 
ideal model, which we have found displays a wide range 
of SDW state, which is entirely non-classical and diffi- 
cult to obtain from spin-waves. The SDW state is very 
naturally related to magnetization plateaus, which have 
been widely discussed, mostly in the spatially isotropic 
regime. This is discussed in detail in Sec. |VII D[ 



C. Open Questions and Parting Remarks 

Our study has resolved most of the main questions re- 
garding the phase diagram of CS2CUCI4. However, there 
are still some smaller details which remain to be under- 
stood. We have not addressed the regime of magnetic 
fields h of order for fields in the triangular plane. 
There, the ground state should involve a non-trivial bal- 
ance between D and inter-chain exchanges. Evidently, 
this gives rise to additional phases for fields along the c 
axis. There is also a discrepancy between the measured 
incommensurability (go) in the high-field cone state for 
this field orientatiorP^ and theoretical expectations. To 
check whether this discrepancy might be explained by 
the additional DM interactions considered in this paper. 
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we evaluated in Appendix[G]the single-magnon spectrum 
in the fully-polarized state, including these interactions. 
We have not found any set of parameters which appear 
consistent with the measured incommensurability, which 
is reduced compared to t he ex pected one. Indeed, the 
theoretical result in Eq. (G20) shows that DM interac- 



tions only increase the incommensurability, making the 
problem worse. This might be a point that warrants more 
extensive experimental investigation. 

While we have focused on the application of our meth- 
ods to CS2CUCI4 (and to a lesser extent, Cs2CuBr4), the 
analysis can be applied to other quasi-ld materials. In- 
deed, it appears that the commensurate AF state de- 
scribed here has been observecP^ in the another trian- 
gular lattice material, CS2C0CI4, in Ref.74 This mate- 
rial is a spin- 1/2 XY-like antiferromagnet forming a spa- 
tially anisotropic triangular lattice. In contrast with the 
naive expectation of incommensurate spiral order along 
the chains, experiments^ finds commensurate antiferro- 
magnetic ordering in the absence of an applied magnetic 
field. This finding matches nicely our description here, 
as in CS2C0CI4 pronounced easy-plane anisotropy (the 
estimat#^ is J^/J^'^ = 0.25, where represents in- 
chain exchange between S'^ components of the nearest 
spins) plays a role quite similar to the external magnetic 
field: it enhances XY spin correlations at the expense of 
Z ones (which, however, remain commensurate with the 
lattice). In fact, the behavior of the compactification ra- 
dius as function of XY anisotrop}'^^ is not very different 



from that in a magnetic field, discussed in Appendix [A] 
27ri?2 ^ i-arccos[JVJ^'^]/7r. Thus the results of Sec. ]V| 
provide a natural theoretical explanation of the observed 
commensurate AF phase in CS2C0CI4. 

To facilitate further application of our methods such 
as this one, we have described them here in sufficient 
detail that they could be readily applied to other prob- 
lems. We hope that the level of detail exposed here 
serves to amplify the tremendous power of the quasi- 
one-dimensional approach, which allows real quantitative 
contact betweem the microscopic spin Hamiltonian and 
universal long-wavelength physics, while including at the 
same time strong ffuctuations and frustration. 
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Appendix A: Notations and parameters 



Throughout the paper the following conventions are used: coupling constants of quantum Hamiltonians written in 
terms of spin densities 5^'^ are denoted as 7. When these terms are expressed in terms of bosonic fields and 0, 
the corresponding couplings change into 7. Coupling constants of various interaction terms of the effective two- and 
three-dimensional Hamiltonians, expressed in terms of the classical phase 1^^ are denoted as g. 

Table [ll] contains a list of notations in which several perturbations are summarized along with their scaling dimen- 
sions. The scaling dimensions are functions of 2ttB? . 



The parameter 27ri?^ as a function of the magnetiza- 
tion M is obtained by solving the Bethe ansatz integral 
equations.l^^K^ Figure [12] shows the numerical data ob- 
tained from P = 2t:R in Ref. [75l Near the saturation 
magnetization, M ~ 1/2, one can solve the integral equa- 
tion analytically and obtain 



M 



(Al) 



In the opposite limit of zero magnetization, M ~ 0, 2ttR^ 
is well fitted by the following function: 



27ri?^ 



1 



1 



21n(Mo/M)' 



(A2) 



where Mq = ^J^/{tic). An abrupt, inverse-log, deviation 
from the SU{2) value is due to the marginally irrele- 



vant current-current interaction typical for the spin- 1/2 
Heisenberg chain.*^ 

The relation between the magnetization M and the 
magnetic field h has been discussed previously in ReffMl 
and we briefiy describe the result here. As discussed 
below (25), the interchain interaction J' increases the 



energy of the system of coupled chains by 2J' M'^ which 
results in suppression of the two-dimensional M{h) curve 
with respect to the one-dimensional Mi{h) curve for a 
single spin chain, at fixed external field h. Using the 



relation h = 



dE 
dM 



, we observe that the field h naturally 



decomposes into a sum of "one-dimensional" field hi = 
~§M ^ single magnetized chain with energy Ei{M) 
and the inter-chain contribution —AJ'M. As a result, we 
arrive at a self-consistent equation where magnetization 
M{h) of the system of coupled chains is approximated by 
that of the single chain, but evaluated at a shifted field 
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Hydrodynamic rep. 


Bosonized rep. 


Scaling dimension 


7sdw(<S^,2;7r-2(5*^^+l,2;7r + 2(5 + ^.C.) 

7sdw — smo 


7sdw C0S[^((/)^,^ - + 

7sdw — J smo 


1 

27ri?2 


-ijconeS^^^.^dxSy^^^^.^ + h.c. 

_ Jf Icy 

/cone — J 1 ^ 


-^corve{dxOy,z + + COs[/3(^^,^ - Oy^l,z)] 
7cone — J /I3P/Z 


2tiR^ + 1 


7± = J /2 


7± cos[/?(ej,,. - ^,,,.+1)] 




— /o /2 


72 cos[/3(6'^,^ - 6*^+2, ;^)] 
^0 — Ai Jo 




biquadratic fluctuation — generated 
coupling 


-7bqCOs[2/3(^y,2 - Oy+l,z)] 




/?c(-l)^Jy% 


dc{-iyd:,Oy,z 


1 



TABLE II: List of notations. 
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FIG. 12: Parameter 27tR^ as a function of the magnetization 
M. Numerical solutions in Ref. |75] sho wn b y the dot s ar e 
compared with asymptotic solutions Eq. (Al) and Eq. (A2). 
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FIG. 13: (Color online) Magnetization M versus magnetic 
field h. Dashed (red) curve shows magnetization of a single 
Heisenberg chain, based on t he d ata in Ref. .75, Solid (blue) 
line shows mean- field result (A3). 



h - AJ'M: 



M{h) = Mi{h- WM) . (A3) 

This equation is easily solved numerically using data of 
Ref. |75 and results in M{h) curve going below Mi{h) 
for all h. An essentially identical result is obtained if one 
uses the interpolating formula Mi{h) = 7r~^ arcsin[l/(l — 
7r/2 + 7r//i)], suggested in ReflU This approximation pre- 
dicts saturation field /igat = 2 J + 2 J' which is very close 
to the exact 2d result = 2 J + 2 J' + {J'f/{2J). 

It is also worth noting that while Mi approaches satura- 
tion with an infinite slope, dMi/dh ~ (2 — /i)~^/^, the 
two-dimensional curve is characterized by the finite slope 
1/(4J0, so that M{h ^ /isat) = 1/2 + (/i - /isat)/(4J0. 



Appendix B: Symmetry analysis 



Here we consider the full symmetries of CS2CUCI4, and 
some of their consequences. Most importantly, we derive 
the possible DM vectors of the on-chain and diagonal 
bonds. The direction of the DM vectors are constrained 
by the space group symmetry in the crystal. In the ideal 
triangular lattice, the DM vectors must be perpendicular 
to the plane and hence is parallel to the a axis. However, 
this is not true in the real crystal due to the lower sym- 
metry. 
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1. Crystal structure and symmetry generators 

CS2CUCI4 has an orthorhombic crystal structure with 
space group Pnma^ The lattice parameters are a = 
9.65A, b = 7.48A, and c = 12.35A at 0.3 K. The unit ceh 



14 



contains four independent Cu^+ ions as shown in Fig. 
The locations of the ions are given hy R-\-S a {(^ = 1, 2, 3, 




^4 

"2W 




FIG. 14: Left: Four independent Cu spins (labeled by 1,2,3, 
and 4) in the unit cell of CS2CUCI4. The spins 1(2) and 3(4) 
lie in the same a-c plane. Right: Layered-triangular lattice 
formed by magnetic sites. Thick, dotted, and broken lines 
indicate J, and J" exchange couplings, respectively. 



and 4), where 



R = la-\- mb + nc 



(Bl) 



denotes the location of the unit cell and 



^1 

^4 



(1 - zo)a + |6+ (1 - yo)c, 



with zq = 0.23 and yo = 0.42. In the simplified nota- 
tion used in Eq. the spins 1,2,3, and 4 correspond to 
^x,y,z, 'S'^+i 5 ^x,y,z+i, and S^^i^y^-^^^^i, respec- 
tively. Note that m in Eq. (Bl) is the coordinate along 
the chains, i.e., x. For fixed /, 1 and 2 spins constitute 
one triangular layer while 3 and 4 another layer. On the 
other hand, 1 and 3 spins lie in the same a-c plane which 
is a mirror plane through the midpoints of the on-chain 
bonds (see Fig. 



14). The same thing holds for 2 and 



4 spins. Therefore, according to Moriya's rule, the DM 
vector on the chain bonds must lie in the a-c plane. In 
contrast to this, there is no symmetry constraint on the 
DM vectors on diagonal bonds. 

We shall next determine the pattern of relative signs 
of the DM vectors, which requires a more elaborate ar- 
gument of symmetry. The Pnma space group has three 
kinds of symmetry operations, A, 5, and reflection apart 



from Br avals lattice translations. We consider the trans- 
formation property of spins under A and B. Let us de- 
note by Sa{l^ TTi^ n) the spin at the position S^x. The 
spins are transformed under A as follows: 



A : 



5'f (/, m, n) ^ S^{-1 - 1, -m - 1, -n - 1) 
S^{1, m, n) ^ S^{-1 - 1, -m - 1, -n - 1) 



(B2) 

where ji = a^b^ c. Since A is the inversion operation as 



shown in Fig. 15, the spins do not change sign. The 
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FIG. 15: Symmetry operations A and B. The locations of 
unit cells are indicated by {l,m,n). 



symmetry operation 5, graphically shown in Fig. [Tsj 
corresponds to a tt rotation about an axis parallel to a. 
Under it, two of three components of the spin change 
sign. Therefore, one obtains 



B : 



Si {I, m, n) ^ 7^6'3 (/, — m, — n) 



82(1^ m, n) ^ JfiS^{l^ —m — 1, —n — 1) 
where 75 = 7c = —1 and 7^ = 1. 



, (B3) 



2. DM vectors 

We can determine t he r elative signs of the DM vectors 
using Eqs. (B2) and (B3). Let us first consider compo- 



nents of the DM interactions on the on-chain bonds: 

4 

^{Da,aa + Da,cc) • Sa{l, m, u) X Sa{l, m + 1, n) (B4) 
Applying A and translations, we find 



^l,c = — ^4,C7 ^2,c = — ^4,c 

while applying B and translations, we find 



(B5) 
(B6) 



(B7) 
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-5'" 



is simply the inversion tranformation A, from Eq. (B2 ) in 
Sec. |B 1[ In the notation of the main text, this operation 
reads 



(B12) 



Clearly, this operation preserves all components of the 
uniform magnetization, and hence is a symmetry for an 
arbitrary applied magnetic field. 

For a magnetic field along a, the B operation given in 
the previous subsection, Eq. (B3), can be chosen as the 



TABLE III: Point group operations in the Prima space group. other symmetry generator. It can be written as 



These six relations determine the relative signs of Da and 
Dc as in Eq. (16). Let us apply the same technique to 



analyze the DM vectors on the diagonal bonds. Since 
each site lies in the mirror plane, we can assume the 
coupling of the following form: 

+ D21 ' S2{l,m,n) X Si{l,m,n 1) 
+ (1^3,2^4) 



with 



D 



± 

a(3 



(B8) 



(B9) 



Here we have used the fact that 5*^ and 5*^ change signs 
while does not under refiection through a mirror 
plane. Applying A and translations, we find 
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D 



43- 



(BIO) 



On the other hand, applying B with translations, we ob- 
tain 



D' 



12, ( 



-D' 



43, a' 



12,6 



^43,6' ^12, ( 



D43,c- (Bll) 



Using Eqs. ([B10| and (|B11|, one can show that relative 
signs are given by Eq. (|17|), which is shown in Fig. [g] 



3. Symmetry with in-plane magnetic field 



(B13) 



with Va = I, Vb = Vc = -1. 

However, operation B does not keep either the b or the 
c-axis magnetization invariant. The second independent 
operation should be chosen differently for these field ori- 
entations. For a magnetic field along 6, it can be taken 
as a refiection in an a-c plane: 



(B14) 



where (a = Cc = -1 and Cb = 1- 

For a field along c, it can be taken as a refiection in an 
a-b plane: 



^ ' ^x,y,z 



(B15) 



where rja = rjb = —1 and rjc = 1. 

Using the bosonization formulae, Eq. ([24|, one can de- 
duce the transformation of the boson field Oy^z{x) under 
these symmetries. One finds, for the inversion operation. 



■ X) 



213 



(BI6) 



For the field along b, the reflection gives 

C : Oy,,ix) ^ Oy,A-x) - ^[1 - (-1)"]- (B17) 

For the field along c, the corresponding refiection instead 
gives 



TT 



(B18) 



To fully determine the spin structures, it is sometimes 
useful to have a more detailed understanding of the resid- 
ual symmetry in a field. The uniform applied field ob- 
viously preserves the translational symmetry of the lat- 
tice, but breaks time reversal symmetry. Its effects upon 
the point group operations are less obvious. In general, 
the Pnma space group contains 7 non-trivial (and one 
trivial) point operations, which are not all independent 
(see Table HI). Of these 7 operations, 3 preserve any one 
component of the magnetization, and of these 3, two are 
independent. Therefore the point group symmetry in the 
presence of the magnetic field is generated by just two op- 
erations, which have Ising character. The first operation 



Appendix C: Breaking of U{1) spin rotation 
symmetry 



In the discussion of Sees. Vl|VI[ the overall phase angle 
of the spins in the plane perpendicular to the magnetic 
field remained arbitrary. This refiects the U{1) symmetry 
of spin rotations about the field axis, which is present in 
the effective Hamiltonian having dropped DM terms with 
DM-vectors perpendicular to the field. While we expect 
this to be an excellent approximation, it is not exact, 
and the weak effects which we have neglected should re- 
move this artificial invariance. This is appropriate for a 
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crystalline system with only discrete symmetries. In this 
Appendix, we use symmetry analysis to determine how 
this symmetry breaking occurs, for different field orienta- 
tions. We will not attempt to determine the microscopic 
origin of these effects here, which might be, for instance, 
symmetric exchange anisotropy, or fluctuation-generated 
interactions. Instead, we ask what terms might arise in 
the energy as a function of the remaining parameters 
describing the orientation of the spins in the plane per- 
pendicular to the field. To do so, we must consider the 
reduced symmetries of the system in the presence of the 
magnetic field. 



1. Cone state 

In Sec. |IV[ we obtained a cone state for arbitrary mag- 
netic fields (below saturation) along the a axis. Consid- 
ering the standard model, this incommensurate state has 
both a U{1) spin-rotational degeneracy and a pseudo- 
U{1) phason degeneracy. The latter degeneracy is pro- 
tected by translational invariance, but the U{1) rota- 
tional symmetry is an artifact, violated for instance by 
the Dc term and other DM interactions with DM- vectors 
along b or c. To study the breaking of the U{1) spin rota- 
tional symmetry, it is useful to consider the combinations 
6± = 6o±0i. The 6_ field describes the phason mode. 
We focus instead of 0+. 

Apart from translation invariance, the residual sym- 
metries in the field along a are A and B. Under these 
operations, we find 



A: 6+ ^ -6+ + 27r, 
5:6+^ -6+. 



(CI) 
(C2) 



One sees from this that the effective potential should be 
an even, 27r-periodic function of 6+. There are no further 
symmetry restrictions on this potential. In the simplest 
situation, e.g. 14ff(0+) = Vcos0+, this potential has 
a unique minimum (9+ = or tt, depending upon the 
sign of V). Thus one expects the artifical U{1) rotational 
symmetry about the a axis to be completely lifted, and 
the only degeneracy of the cone state to be that associ- 
ated with the phason mode. 



2. Antiferromagnetic (AF) state 



a. field along b axis 

For a field along the b axis, the remaining symmetries 
are the inversion and reflection operations, denoted A 
and C given in Eqs. ( |B12|B15D , in Appendix |B3j Using 
the formula in the appendix, and Eq. (80), we find that 
under these operations. 



A: a ^ a, 
C : a ^ -cr. 



e^6+-(l + a), (C3) 



9^6. 



(C4) 



In addition to these symmetries, under translations by 
one unit along x, one has 



6^6 



(C5) 



From this, we may construct the simplest energy func- 
tion which depends upon O and a, in the spirit of Lan- 
dau theory. From C, Eq. (C4), we see that it must be 
independent of a. From A and T, Eqs. ( C3|C5 ), it must 
be periodic in O with period tt. Importantly, there are 
no other symmetry restrictions. Hence the most general 
energy is of the form 



v{e) 



oo 



. cos(2nO — an)- 



(C6) 



Generically, since there are no restrictions on a^, such a 
potential V(0) has two inequivalent minima (e.g. if we 
take the simplest form with = for n > 1, located 
at O = ai/2,ai/2 + tt) located at points of no partic- 
ular symmetry. This means that the transverse (to the 
magnetic field) components of spins do not lie parallel 
to either the a or c axes. In total, one therefore expects 
4 ground states, with (a, 6) = (±1, 9o), (±1, 6o + tt). 
Physically, the four states are obtained from one another 
by the arbitrary choice of global sign for the transverse 
components of the spins on the even and odd chains, sep- 
arately. 



b. field along c axis 

In this field orientation, the good symm etry operations 
are A, T and given in Eqs. ( |B12|B15D of Ap pendix [B] 
A and T continue to act as in Eqs. ( C3|C5 ), while D 
leads instead to 



The breaking of the spin-rotation symmetry about the 
field axis is crucial for determining the precise nature of 
the spin structure in the AF phase. Since this phase 
is commensurate, and the Hamiltonian has only discrete 
symmetries, we expect only a discrete ground state de- 
generacy. This will fix the angles of the spin projections 
into the plane transverse to the magnetic field. We want 
to construct an effective potential which depends upon 
the parameters of the AF state, O and a. 



D : a ^ -cr, 



6 ^ e + TT. 



(C7) 



Combining the D and T operations, one sees that the 
energy must be independent of a, and then using A or 
one obtains again the effective potential for O in form 
of Eq. (C6). Thus, the ground state degeneracy (four- 
fold) is the same as in the AF state for fields along 6, 
and the transverse components of the spins do not point 
along the a or 6 axes. 
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3. Incommensurate Phase 

Here we must reconsider the vahd symmetry opera- 
tions - A, T, and D - when acting upon the spin struc- 



ture in Eq. (100). This structure is parametrized by two 
angles, 6o and 6i, which describe the spins in the even 
and odd chains, respectively. Equivalently, we can choose 
the combinations Q± = 0o ± Oi. We find 



A:e+^e+ + 27r, 6_^-6_ 
D:6+^6++7r, 6_ ^ e_ 
T : e+ ^ 6+ + 27r, 6_ ^ e_ 



(C8) 
(C9) 
■2qo. (CIO) 



Note that the transformations of the B_ field under A 
and T involve shifts by multiples of the incommensurate 
wavevector qq. Under multiple actions of such shifts, any 
value of 6_ can be approached arbitrarily closely (due to 
the 27r periodicity) . Thus there is no potential which can 
pin the values of 0_. This is the reason for the gapless 
"phason" mode in the IC phase. 

By contrast, the effective Hamiltonian can certainly 
depend upon 9+, refiecting the lack of rotational invari- 
ance about the c axis. In general, due to the action of D 
(which is most constraining), the effective potential must 
be a TT-periodic function of 6+. There are, however, no 
other constraints. Since shifts of either 6o or 6i by 27r 
have no physical significance, 0+ is itself defined only up 
to 27r. This implies that there should be two discrete sets 
of IC solutions, described by 6+ = 6o, Oo + tt, and 6o 
should generically take an incommensurate value with no 
special symmetry. 



Appendix D: Chain mean field theory 

1. SDW order 

We start by applying CMFT to the ideal 2d model 
Eq. (27) and consider first SDW order at finite tempera- 



ture. Thus we write 



^i,sdw J ^^7sdwCos[27r(^^„ 



>,+i,,)//3]. (Dl) 



It is convenient to shift (j) fields slightly, <py^z ^ ^^y^z + 
(—1)^/3/4, so as to change the sign of the interaction in 
the equation above. In terms of the shifted fields the 
minimum corresponds to a uniform configuration (j^y^z = 
(j)Q{z) for each since the layers are decoupled. It is clear 
that the inter-chain interaction can next be written as 
(^y,z • (^y+i,z where tr^,^ = (cos[27r^^,^//3], sin[27r0^,^//3]) 
is the SDW order parameter vector describing chain 
{y^z). Chain mean-field consists in a self-consistent as- 
sumption that SDW order spontaneously develops at 
some critical temperature Tgdw, below which the SDW 
order parameter acquires a finite value along some arbi- 
trary direction in the SDW plane. For concreteness we 
choose this direction to be along x-axis, {(Ty^z) = (V^^O)- 



This choice corresponds to = in Eq. (28). With 
these approximations we have 



-2tjj%dw ^ dx cos I 
y,z 



1 



(D2) 

where the factor of 2 arises from the coordination number 
of chain y in the layer z. The self-consistent condition 
reads 



^ = (C0s[^^2^,^])sdw, 



(D3) 



where angular brackets denote finite-temperature aver- 
age with the sine-Gordon Hamiltonian of the single 
chain: 

H^l = Jdx'^{{d,ey,^)' + {d,4>y,z)')+H^L^. (D4) 



The right-hand side of Eq. (D3) is evaluated perturba- 



tively in powers of vanishing and the leading order 
result is 



ijj 2V;7sdwXA,dw(<? = 0,cJn = 0;T). 



(D5) 



Here we defined the momentum and frequency depen- 
dent susceptibility^ XAiQ^^n^T)^ of the vertex operator 
Oa = cos(V47rA(/)) (or Oa = cos(V^7rA^), which gives 
identical results) at temperature T, in the canonical free 
boson theory, Eq. 

/OO nl/T 
-OO -'O 

x(Oa(x,t)Oa(0,0))o.(D6) 

The subscript reminds us that that this is evaluated 
in the free theory. This susceptibility, in various limits, 
plays a central role in the determination of critical tem- 
peratures within CMFT. It is evaluated at the end of this 



appendix. Sec. D 7 Here, we need the SDW susceptibil- 
ity, for which Ag^w = l/47ri?^. At the critical temper- 



ature Eq. ( |D5[ ) acquires a non-trivial solution, ip 0, 

(D7) 



resulting in the implicit equation for Tgdw* 

1 = 27sdwXAsdw(<? = 0,uJn = 0;Tsdw)- 
Using Eq. (D55) from Sec. |D 7 a] we obtain 



^ 27rTsdw ^ 



2-2A,, 



r(i-A,dw)r^(Asdw/2) 
^^^-r(A,dw)n(i-A,dw/2) ^""'^ 



1 + Asdwr(Asdw - l/2)/(v^r(Asdw)(l - Asdw)) 



where Agdw = 7r%dw/v = 7rAlsm{S)jyv. This result, 
Tgdw as a function of magnetization M, is plotted in Fig. 
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2. Cone order 

Here we consider the cone (twist) ordering instability 
of the ideal 2d model Eq. (27). The cone Hamiltonian is 
given by 



^l,c 



-7c 



X cos 



(D9) 



where^cone = -^'^3/^/2. The spatial derivatives in 
Eq. (D9) require a generalization of the procedure de- 



scribed in subsection |D 1[ We begin by shifting the 
fields by a position-dependent phase corresponding to a 
wavevector shift g'o, the magnitude of which is to be de- 
termined later self-consistently. Thus 



This shift transforms Eq. ( |D9| ) into 

^<707cone 



(DIO) 



Hi, 



5/ 



(Dll) 

where we have neglected as subleading dxO terms. The 
transformed Hamiltonian Eq. (Dll) is now of the form 



Eq. (|D1| and c an be manipulated similarly. However, the 
shift Eq. ( |D1Q ) has modified the free boson Hamiltonian 
Eq. ^ into 

(D12) 



Introducing the order parameter 

= (cOs[/3l9^,^])cone, 

where the average is over the Hamiltonian 



(D13) 



H^^^ = Ho- 2^7eo„e Y.J^^ cos[0y^,] , (D14) 



and 7cone = 2(2o7cone//3 = qo J'Al. As before, expanding 
the expectation value in Eq. (D13) to leading order in 7, 



and assuming ip 0, gives the condition for the critical 
temperature 



l = 27coneXA±((/ = 0,0;Te 



(D15) 



where tilde on the susceptibility indicates that it is to 
be calculated using the free but shifted Hamiltonian 
Eq. (D12). However we can now undo the shift Eq. (DIO) 



and transform back to the original fields. As a result, 
one obtains the identity 

XA± {q = 0, 0; T,) = xa± (^o, 0; TJ. (D16) 
Then Eq. ( D15[ ) becomes 

l=27coneXA±(^0,0;Te). (D17) 



This wavevector dependent static susceptibility is well- 
known, and given in Eq. (D51). Using it, Eq. (D17) 



can be solved by maximizing the right-hand-side with 
respect to qo at given T, and then finding the maximum 
Tc = Tcone (this is equivalent to choosing the wavevector 
qo for which Tc is maximum). Expressing all quantities 
in terms of dimensionless variables, r = vqo/{4:7TTcone) 
and s = 27rTcone/'^7 we obtain a system of two equations 



4 Im^(^ + ir) 



a-2A± 



27r sinh[27rr] 



cosh[27rr] — cos[7rA±] 



27coner(l- A±) 

7r(3v r(A±) 



r(- 



x(cosh[27rr] - cos[7rA±]). (D18) 

The resulting Tcone is plotted in Fig. (T] 

It is worth mentioning here that the outlined calcula- 
tion can be done by keeping track of lattice as well, so 
that spatial derivatives in Eq. ( D9 ) become lattice differ- 
ences. Following this route (see for example Ref. 18) one 
again arrives at Eq. ( D17| ) but with the coupling constant 
given by 7cone,iattice = 2J'Alsm[qo/2]. This difference, 
qo 2sin[g'o/2], does not affect the outcome as in the 
regime where CMFT is applicable the ordering vector re- 
mains small, qo ~ J' / J <C 1. See section D6 for more 
discussion. 



3. Inter-layer interaction J'' 

Here we consider the fate of SDW and cone orders 
in the presence of inter-layer coupling J" . As discussed 
in Section |IV[ the inter-layer interaction is a strongly 
relevant perturbation which should be accounted for in 
CMFT. 



a. SDW order 



Consider the SDW channel first. Eq. (Dl) should now 
be complimented by 



^2,sdw — ^ / dx (5^,2;7r-2(5*5^,2+l;7r+25 

y,z 

+*^y,2;7r+25*^^,2+l;7r-2(5)- 



(D19) 



This is simply a rewriting of the 7^' term in Eq. (57). Its 
bosonized form is 



^2,sdw = l^WYl J C0S[2^((/)^,, 



(D20) 

The total SDW Hamiltonian is obtained by adding 
Eq. (Dl) and Eq. (D20). Both terms can be made neg- 



ative by a shift 0^,^ ^ 0^,^ + ((-1)^ + (-l)'')/^/^. Fol- 
lowing the same steps as in Sec. |D 1[ we find that Tc is 
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determined by an equation of the same form as Eq. {D7), 
but with the replacement 



7sdw ^ r,dw = 7sdw + A\J"I2 = A\{J' sin (5 + J" jl). 

(D21) 



Fig. [16| shows that Tgdw is mildly enhanced by J" at low 
magnetization M. 
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FIG. 16: (Color online) SDW ordering temperature Tsdw 
(dotted (red) line) of the ideal 2d model, as obtained from 
Eq. (D7). Solid (blue) line: the same but with inter- layer J" 
accounted for, see Eq. (D21). 



h. transverse/ cone order 

Interlayer exchange J" strongly enhances the trans- 
verse (to the field) order (underlying the cone, AF, and 
IC states) for magnetic fields along the h and c axes, for 
which the DM interaction D is ineffective. In this case, 
we need to account for the (7^^) term represented in 
Eq. ([71]), which is bosonized in Eq. ([72|). To bring the 



latter into canonical form for CMFT, we change its sign 
by a simple shift dy^z ^y.z^^^l which does not affect 
any of the manipulations in subsection D 2 Transforming 



next to dy^z as in Eq. (DIO), we arrive at the modification 



of Eq. (D17) where 7, 



is replaced by 



J") 



(D22) 



As in Sec. |D 2[ we obtain two equations for and T by 
maximizing the right-hand side of the modified Eq. (D17) 



with respect to q^^ and using the equation itself. The re- 
sult rea ds, in terms of dimensionless pair (r, s) introduced 
in Eqs. (IdT8|), 



4 Im^(- 



2 

.2-2A± 



ir) 



27r sinh[27rr] 



2J^ 



cosh[27rr 

Al r(i 



— cos[7rA± 



2J'sr + J'^ 



2tiv r(A±) 
x(cosh[27rr] - 



cos[7rA± 



2J'sr + J" ' 

4 



• ir) 



(D23) 



The cone ordering temperature Tcone is plotted in Fig. [17 
We observe that the interlayer coupling enhances Tcone 
dramatically, and even leads to a substantial Tc when 
J' = 0. The reason for this is simply that the non- 
frustrated nature of the interlayer exchange leads to an 
appreciable inter-chain coupling of transverse spin com- 
ponents even when q^ <C 1, as Eq. (D22) shows. Indeed, 



in the AF and IC phases, we also have transverse order- 
ing, but J' does not actually contribute to and the 
plot with of Tc{J' = 0) is relevant in those cases. 
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FIG. 17: (Color online) Critical temperature of cone order- 
ing as obtained in the ideal 2d model, Eq. (D9), (solid (blue 
line) and in the 3d model (dashed (green 
Eq. (|D23|). The dot-dashed 



line) as given by 
magneta) curve represents the 
critical temperature for Tn of the interlayer-driven order, rel- 
evant for the AF an d IC phases, which we obtain by setting 
J' = in Eq. (|D22). 



The effect of interlayer exchange J" on the two orders 
can now be compared. Notably, despite its smallness - 
J" = 0.045 J - interlayer coupling completely eliminates 
SDW order in a system of weakly coupled layers, T^dw < 
Tcone foi" all magnetizations from to 1/2, see Fig. 18 



4. Field along a axis 

Here we describe how to calculate Tc in CMFT in the 
presence of the DM interaction D for a field along the a 
axis, which is the arrangement considered in Sec. IV of 



the main text. Here J'' is unimportant, as we will see. 
Bosonization of the DM term in Eq. (57) gives 



^2, dm 



/ 



dx {-iy2DAl cos[(3{0y^z - Oy+i^, 



(D24) 

while the cone term is given by E g. (|D9| ), with 7cone = 
J' A\l3/2. We observe (c.f. Section [IV B| where the corre- 
sponding T = state is discussed) that the two interac- 
tions can enhance each other if the sign of the ordering 
vector is correlated with the sign of the DM vector. We 
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FIG. 18: (Color online) Tsdw (solid (blue) line) and Tcone 
(dashed (green) line) for J" — 0.045J. Compare this with 
ideal 2d situation in Fig, [7| where the two orders compete 
strongly at intermediate values of M. 



FIG. 19: (Color online) Ordering temperature with the field 
along the DM (a axis) axis. Dashed (cyan ) curv e: Td{J' — 0) 
due to DM interaction only {J' = in Eq. ( |D27| )). Solid (blue) 
curve: Td obtained with both D and J' interactions present. 



therefore make the layer-dependent shift, which corre- 
sponds to Eq. (p5| , 



ey,,{x) = i-lYqox/p + 0y,,ix). (D25) 

This should be contrasted with Eq. ( |D1Q ), which de- 
scribes the situation without any DM vector. The trans- 
formation in Eq. (D25) makes the competition between 



the staggered DM and interlayer interactions (discussed 
in Sec. IV) obvious, since the argument of the inter- 



layer cosine term, in Eq. (72), acquires a position- 
dependent phase 



= {-iy2qox/p^i 



(D26) 



The resulting oscillations eliminates the contribution of 
J'' to the energy within CMFT. Proceeding as described 
in subsection |D 2[ we again obtain an equation for the 
critical temperature in the same form as Eq. (D15), but 
with 7c, 



replaced by F^^ 



2D). 



(D27) 



Note the great similarity of the above coupling with that 
in Eq. (D22): the two situations are related by exchang- 
ing J'^ ^ 2D. Hence the critical temperature, T^, 
for transverse (cone) type ordering follows from solving 
Eq. (|D23|) with J" replaced by 2D. The result is plotted 



in Fig. |19[ which compares the case of DM interaction 
only Td{J' = 0) with that of general D 7^ 0, J' / 
situation. One observes that J' leads to only a modest 
enhancement of Tjj relative to the case with DM interac- 
tion present only. Note that Figj8]in the main text shows 
the solution with D ^0,J' = 0. 



5. CMFT at T = 



Here we outline calculation leading to Eq. (76). We a\M) 
start by changing the sign of Eq. (72) via a shift: Oy^z 



Oy,z + {-iy7r/{2P). Within CMFT, is replaced by a 
single-chain sine-Gordon Hamiltonian, 



K.0 = -H'J 



dx cos j30{x). 



(D28) 



Here il) = {cos f3 6) is to be determined self-consistently. 
The brackets stand for averaging with the sine- Gordon 
action which, upon rescaling of the temporal coordinate 
r = y/v, reads 

^sG= I dxdy(^{djf ^]^{dyef -2^cospe). (D29) 



Here /i = \IjJ"AI/v. The exact solution of Eq. ( |D29[ ) 
from Ref. [78l gives the ground state energy density 



F(/i) 



1 



rMnan 



(D30) 



which is expressed in terms of the parameter ^ and mass 
gap M. These are determined by 



/32 



1 - /3'2 Stt - ' 



r(/?'2) 



M 



v^r((i + e)/2) 



2r(e/2) 



2-2/3'^ 



7rr(l - /3'2) 

with /3' = l3 /\/8n. Using the obvious relation 

1 dFin) 



(cos/36') 



2 di^ ' 



we obtain, after some algebra, 

'J"Al \ 2-2"il2 



= a'{M) 



tan[7rC/2] 
27r(l - j3'^) 



r(i) 



7rr(l - 



r(/3'2) 



,(D31) 



(D32) 



(D33) 



1/(1-/3") 



35 



The order parameter in Eq. ( 76 ) is related to the self- 
consistent ly calculated ip very simply, ip = As^p. Hence 
the prefactor in (PtgI) follows as 



Cr(M) 



ll/(l-A±) 



(D34) 



The resulting order parameter, ip{M)^ is plotted in 



Fig. 20 as a function of magnetization. Because the ex- 
ponent TiR^ /{2 - 2tiR^) in Eq. ( |D33[ ) is always smah (it 
varies from 1/2 at M = to 1/6 at M = 1/2), one ob- 
serves that ipiyM) is a slow- varying function of magnetiza- 
tion M. The overall non-monotonic shape of ip{M) is in 
excellent agreement with the experimental data reported 
in Fig. 3b of Ref. [161 the order parameter first rises 
with increasing magnetization, reflecting the increasing 
relevance of transverse spin correlations, and then falls 
down rapidly as M ^ 1 /2 as a result of diminishing den- 
sity of magnons. However, the CMFT result shown in 
the figure does suffer one problem, discussed in the next 
subsection, leading it to break down close to saturation 
(where it incorrectly predicts that ip remains finite as 
M 1/2). 
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FIG. 20: The T = order parameter versus M, calculated 
within CMFT. 



6. Limitations of CMFT 

The results of CMFT for the critical temperature and 
order parameter, discussed in the previous subsections, 
exhibit unphysical behavior on approaching the satura- 
tion magnetization: in these calculations, the Tc for the 
cone state remains finite in this limit, as does the zero 
temperature order parameter ip. These features are both 
clearly incorrect, as both Tc and ip must decrease to zero 
as the spins become fully polarized. We will investi- 
gate the breakdown of CMFT in this subsection in more 
detail, and determine the proper scaling theory for the 
vicinity of magnetization saturation. 

First, let us observe the failur e of C MFT more directly. 
Consider the solution of Eqs. (D18) in the A± 1/4 



limit in more detail. First, note that the first equation 
there is readily solved by r = 0(1) which is not par- 
ticularly sensitive to the value of M. This immediately 
tells us that s - {Ap'/v)^ for A± = 1/4. Near the 
saturations^ ^ (^-M)^/^ while v ^ {\-M)J, which 
implies s - {J'/Jf /{\ - M). Hence vs = {J'f/J 
approaches a constant value while the ordering momen- 
tum = 2rs - {J' /Jf/{\ - M) diverges as M ^ ^. 
The latter divergence is a clear indication of the failure 
of CMFT. It can be traced the fact that CMFT is by 
construction an expansion about the Id chain limit. The 
natural parameter of this expansion is J' /v which is sup- 
posed to be small everywhere. This assumption clearly 
breaks down near saturation, where the spin velocity v 
vanishes and the expansion is not justified anymore. 

Physically, near saturation, one has a dilute gas of 
spin flips, which can be thought of as hard-core bosons 
or, in one dimension, equivalently, as spinless fermions. 
Their density n = \ — M determines Fermi-momentum 
kp = 7rn, which in turn determines the (Fermi) veloc- 
ity as V = kF/m. Since the mass m is of the order of 
inverse chain exchange constant J, we obtain the scal- 
ing quoted above, v ^ — M)J . The velocity vanishes 
because precisely at saturation, the hard-core magnons 
possess a quadratic dispersion relation, which is beyond 
the Luttinger liquid paradigm of linearly dispersing col- 
lective excitations. (The scaling of A3 ~ v}^^ follows 
from the fact that the scaling dimension of the spin-flip 
operator is iiR^ = 1/4.) 

To solve the problem correctly we need to start with 
a 2d description, which is actually simple near the satu- 
ration. We define the "order parameter" field ^^(x) ~ 
Sy.^{x), which is just the annihilation operator for a spin 
flip. Because the flips are dilute, we expand their ki- 
netic energy near the bottom of their Id band, and the 
Hamiltonian of the layer is given by 



'2m 



/i)^^(x) 



(D35) 



where the chemical potential ja ~ /igat — h describes the 
deviation from the saturating magnetic field h^at and dots 
stand for the interaction terms (which are irrelevant at 
the 2d critical point). The scaling behavior for small 
J' / J and small deviations from saturation, can already 
be read off from Eq. ( D35[ ), which should be considered 
for this purpose as a (T+l) -dimensional quantum field 
theory. The magnetization density, relative to saturation 
is ^ — M ^ ^"f"^ which scales like an inverse length. Fur- 
thermore, the J' term above scales in the same way, since 
it is missing one x derivative. Thus physical quantities 
should be scaling functions of E! = [J' / J)/{^ — M), as 
claimed in Eq. (30). The overall scaling behavior of Tc is 



determined by the dynamical critical exponent, which is 
z = 2 due to the quadratically-dispersing magnons. Since 
Tc therefore scales as the square of an inverse length, the 
result in Eq. (31) follows. 
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We may understand the result more physically, and in 
particular the b ehavio r of the scaling function J-'(S) by 
considering Eq. (D35) and its consequences in more de- 
tail. First, the behavior for small S corresponds just 
to the CMFT result in the limit of M 1/2. The 
limit of large S is more interesting, and represents 2+1- 
dimensional physics. Transforming to momentum space, 
we find that the magnons have the following dispersion: 



J kr, COS k 



2m 

{kx — m J' cos ky)'^ m{J') 



(D36) 



i\2 



cos^ ky — jl. 



2m 

We expand around the minimum, writing kx = mJ' ^qx^ 
ky = qy, which gives the continuum dispersion relation. 



2 ^ 2m 



(D37) 



The continuum theory describes a free Bose gas with 
anisotropic effective mass. It of course forms a Bose 
condensate at low temperature, which is described in 
the usual way by taking the continuum limit ^y{x) 



AmJ'x 



^(x, with 



iO{x,y) 



(D38) 



Here is the 2d superfluid phase, and Ug ^ ^ — M. For 
the above anisotropic Bose condensate, standard manip- 
ulations give the phase-only effective Hamiltonian (su- 
perfluid kinetic energy). 



^jdxdy [px{dxef + Py{dyOfy (D39) 

Eq. ( |D39 ) describes an anisotropic 2d XY model with 
stiffnesses px = n^/m ~ UgJ and py ^ m{J')'^ns = 
ns{J')'^ / J. Its critical temperature is determined by the 
geometric mean of the two stiffnesses: Tc = 7T^pxpy/2 ~ 
Us J'. This argument describes the large S limit of the 
scaling function ^(S). Note also that this argument im- 
plies that the ordering momentum qo ~ T^one/"^ remains 
small even near the saturation. This we saw already 
in spin- wave theory in Appendix [Gj where we observed 
^0 J'/J- 

Away from saturation, so that M — ^ ^ J' / J ^ the 
CMFT is applicable, and Eq. ( D18[ ) predicts small go 
again. For this reason we expect that q^ is uniformly 
small for all magnetizations, and lattice effects, of the 
kind mentioned in the end of Section pD 2[ remain unim- 
portant in all the regimes considered. For this reason we 
chose to keep only the leading terms of the small wavevec- 
tor expansion (that is, approximate lattice differences by 
spatial derivatives) throughout the main text. 



7. Calculation of susceptibilities 

Here we present some technical details of the evalua- 
tion of susceptibilites used in this appendix. We define 



for convenience the susceptibility in space and time, 

XA(^,r,T) = ((9A(:r,r)(9A(0,0))o.:r, (D40) 

evaluated at temperatur e T. It is straightforward to per- 
form the average in Eq. ( D40 ) with the free boson Hamil- 
tonian Eq. ([7|), see for example, Ref. 79: 

XA(^,r;T) = iexp{-4^A x /(y,r)}, (D41) 



dq 
2^' 



-a\q\ ■ 



Here y = x/v^ a = ao/v is the short-distance cutoff 
needed to regularize the integral, cj^ = 27rTn is the stan- 
dard bosonic Matsubara frequency, and 0<r<l/Tis 
the Matsubara time. The frequency summation is stan- 
dard (we use GR 1.445.2 in Ref. 80) and leads to 



dq 



-aq 



iqy-qr 



iqy-qr 



q 1 



^iqy+qr 



-iqy+qr- 



(D42) 



Next we expand the denominator in series which is eval- 
uated term by term with the help of the identity 



d_qf_ 

q ^ 



-Aq 



-Bq^ 



(D43) 



In this way we arrive at 



47r/ In ■ 



m=l 



-In 



(^In 

(m + aT)2 



{m^aTf 



T\T-iyf 



{m^aTY 
-T\T^iyY 



(m + aTY 



(D44) 



For later use we note here that small (r, y) ~ behavior 
is described by the first term in the right-hand side of 
the above equation. It is easy to check that small (t, y) 
behavior, where t = 1/T — r, is described by a similar 
ln[(a + tY + y'^] term (which is contained in m = 1 con- 
tributions). 

Focusing for the moment on the regime where \t\ ^ 
a (i.e. including the small r limit but not the limit 
near r = 1/T), we next observe that m ^ contain 
no singular dependence on a, which allows us to set 
a = there. The well-known identity ln[sin(x)/x] = 
ln[l - x^/in^k^)] leads us to 
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{a + r)^ + y'^ sin[7rT(r - iy)] sin[7rT(r + iy)] 



7r2T2(r2 + y2 



Hence we obtain, when \t\ a, 



(D45) 



2xa{x, r]T) = I 



-r 



{a + r)2 + y^ 



' sm[7TT{r -iy)]sm[7TT{r ^iy)]i ' 
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which is imaginary time version of the expression ob- 
tained in Ref. (STl The appearance of the cut-off a in the 
temporal r direction is a generic feature of the bosoniza- 
tion techniqu e (see Ref. 82 for more examples). The first 
factor in Eq. ( D46 ) clearly carries the information on the 
small r limit. To account for both r ~ and r ^ 1/T 
properly, we must include another similar factor, 



2xa{x, t;T) = { — 
I- (a 



T2 + y2 



.2 ,,y2 



X 



, . (D47) 

sin[7rT(r - iy)] sin[7rT(r + iy)] J ^ ^ 

Eq. ( D47[ ) then is correct for the full imaginary time in- 
terval. To separate the short and long time behaviors of 
Xa(^,t;T) we approximate it as 

Xa{x, r; T) = xi(x, r; T) + xli^. r; T). (D48) 

Here the first term describes long-distance behavior. 



X 



sin[7rT(r — iy)] sin[7rT(r + iy)] J 



xi(x,r;T) = { 

(D49) 

which gives the naive limit of Eq. ( |D47| ) when a ^ 0, 
valid away from the endpoints. The second. 



(D50) 



accounts for an important difference in the short-distance 
behavior of the full Eq. (D47) and the approximate 



Eq. (D49) expressions for Xa(^, T). Observe that the 
short-distance behavior is not sensitive to temperature 
as it takes place on the scale determined by a while the 
long-distance one describes correlations on a much-longer 
thermal scale 1/T. 

For CMFT, we require certain limits of the Fourier 



transform of the susceptibility, Eq. (D6), or 



XA{q^^n;T) 



(D51) 
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a. SDW case 

For the SDW case, we need the static (cj^ = 0), zero 
momentum {q = 0) limit of Eq. ( D51[ ) in the range 
1/2 < A < 1. Here we must take some care to keep 
the short and long time contributions separate. We split 
Xa(0,0;T) = Xa(0.0;^) + X1(0,0;T), Fourier trans- 
forming separately Eq. ( |D49 ) and Eq. ( D5Q ). The first 
contribution is rather standard and given by 



xi(0,0;T) = 



(D52) 



TTVa 



(27rTaj 



2A-2r(i- A)r2(A/2) 
r(A)r2(i- A/2)* 



Observe that Eq. (D52) diverges as 1/(1 — A) when 
A ^ 1. In the case of the SDW order this limit cor- 
responds to the behavior near saturation where A ap- 
proaches 1 at the full magnetization, M = 1/2 (see 
Table. [l| and ( |A1| )). The divergence is not physical 
and stems from the incorrect short-distance behavior 
of Eq. ( |D49| ). It is compensated by Xa(0'05^)' which 
yields t wo id entical contributions from the two terms 
in Eq. (|D50|). Substituting r = at^y = az and using 



dz/ilXz^)^ = v^r(A - l/2)/(2r(A)) we arrive at 



X1(0,0;T) = 
va^^/^T{A-l/2) 
(2 - 2A)r(A) 



((1 + ^) 



2-2A 



L 



2-2A 



(D53) 
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where L = l/{aT) 1. Simple calculation shows that 
this expression diverges logarithmically near A = 1/2, 
XAsdi/2(^' ^' ^ In(Le), which however represents a 



small subleading correction to Eq. (D52), which diverges 



linearly, XARdi/2(^' ^' ^ ^' ^^^^ region. Thus 
XAf^i/2(^' ^' safely neglected. Near A = 1 
limit things are different: here Xa^i{^^^'^T) results in a 
large T-independent contribution 

XA«ilu,u,Jj (2-2A)r(A) ~ 2(1 -A)' 

(D54) 

Similar short-distance correction can be found in 
Refs. I81|83l Collecting both contributions we finally ob- 
tain, for 1/2 < A < 1, 



Xa(0,0; T) 



■{( 



27rTa 1 



. 2A-2 



r(i- A)r2(A/2) 



2 IV J r(A)r2(i- A/2) 

r(A-i/2) ^ 
v^(i- A)r(A)/- 



(D55) 
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At A = 1 this expression reduces to 

Xa=i(0,0;T) = -TTva^ ln[27rTa], 

which is free of unphysical (1 — A)~^ divergence. The 
resulting Tc(A = 1) is exponentially small in v/%dw 
ratio.EHl 

h. Cone order 

For the CMFT treatment of the cone state, one re- 
quires the static susceptibility at non-zero wavevector, 
with A = A± = 7tR^. This is always less than or equal 
to 1/2, making the short-time corrections in Eq. (D50) 



negligible. Therefore we may directly Fourier transform 
only the long time term, Eq. (D49). The result is well- 
known (see, e.g. Ref. [18]): 



XA(g,0;0) 



27TTa 



x2A-2r(l- A) 



r(A) 



r{A/2Xivq/{47rT)) 



r{l- A/2Xivq/{47rT)) 



(D57) 
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Appendix E: Generation of biquadratic interaction Here we introduced the shorthand = e^^^^y . To first 



In this appendix, we detail the generation of the bi- 
quadratic interaction, Eq. (83 ). We use a standard Wilso- 
nian RG, in which one derives the low-energy theory by 
integrating out high-energy modes. 

We begin by passing from the Hamiltonian formulation 
to the (Euclidean) Lagrangian one, integrating out the 
conjugate field in the path integral. Furthermore, we 
rescale the temporal direction, introducing xq = vr, Xi = 
X, in order to render the free action rotationally invariant 
in the x = (xq^xi) plane. The free action of the fields, 
corresponding to Eq. ([7|), is 

1 



^0 



1 r' 



Oy{k)Oy{-k). 



(El) 



Here we have introduced a momentum space cut-off A, 
which is actually of 0(1) (the lattice spacing). Further- 
more, throughout this appendix, we have suppressed the 
z index of the chains to simplify the formulae. Note that, 
due to the change from r to xq, the interaction terms be- 
come perturbations with dimensionless couplings, given 
by original ones divided by v. 

The RG proceeds in the standard way, by progressively 
integrating out modes within a shell of width dA = Adi 
near the cut-off, thereby reducing the latter to a "run- 
ning" cut-off A£ = Ae~^, which defines the logarithmic 
scaling variable i G {0,cx)}. Following the convention 
used in the main text, we do not perform any iterative 
rescaling of length and time scales, thereby allowing the 
cutoff to "run" to increasingly smaller value as the RG 
proceeds. 

Formally, the integration of modes is accomplished, in 
each iteration, by writing 

(E2) 



9< 

^2/ 



9> 



where the "slow" field Oy contains non-zero Fourier com- 
ponents with k < Aij^di = Aie~^^ 



and the "fast" field l9> 



contains the remaining ones with A^j^di < k < A^. We in- 
tegrate out the fast field at each iteration, perturbatively 
in the interactions. After this, we relabel 6y Oy^ which 
then defines the theory at the reduced cut-off A^j^dl- 

At zeroth order in the interactions, the free action 
renormalizes trivially, since the slow and fa st m odes are 
decoupled. It remains in the form of Eq. (El), with A 
replaced by A^. To first order, the perturbations ^1,2,3 
renormalize very simply, according to their scaling di- 
mensions. We illustrate this explicitly for the cone/twist 
interaction under consideration here. Ignoring the SDW 
term, we have the action corresponding to H[ in Eq. (27): 



S'l = -7cone ^ j dxdr {d^Oy + d^Oy^i) cos I3{0y - Oy^i 



(E3) 



order in the RG, we use Eq. (E2) and average Eq. (E3) 
over the fast fields using the free action. Since the fields 
are decoupled at each the two factors average in- 
dependently. One has 



exp[ 



> 

2 pAi 



d^k 1 



(E4) 

where A± = /3^/(47r) = ttR^ is just the scaling di- 
mension of the field. Le tting 0^ 0, we see that 
Eq. ( |E4[ ) applied to Eq. ( |E3| ) simply multiplies 7cone by 
the constant e~^^±^^. Hence, we have 7cone(^ + d£) = 
(l-2A±d^)7cone(^),or 



-2A±7, 



This of course integrates to 

7cone('^) — 7cone(0)6 



.(0) 



2A± 



(E5) 



(E6) 



The same treatment holds for the interlayer couphng 



^'l (see Eq. (72): 



7±W 



2A± 



(E7) 



As we have discussed, this is the most strongly relevant 
interaction for fields in the b-c plane, which we consider 
here. The RG can be considered perturbative provided 
the dimensionless rescaled coupling, j^J_{£)/v^ remains 
small compared with the typical value of the bare action 
at the corresponding scale, A|. This fixes the value of the 
cut-off, A^', at which the coupled chains form correlated 
a-b planes: 



^(A7A)2^± 

V 



: (A^^ 

Solving for A'^, and using 7^ = A^J'^ 



A' 



^jMiy/(2-2A±) 



(E8) 



(E9) 



where we have used that the bare cutoff A is 0(1). 
Thus the corresponding spatial scale which determines 
the renormalization of all interactions is given by ^" = 
1/A'' - (v/J'')i/(2-2^±), in agreement with Eq. ([74). 

Let us now, finally, generate the biquadratic term. 
This occurs as a second order contribution of 7cone to 
the effective action. Expanding the action in powers of 
this term, we get in second order {Z = J e~^^[l + S'(i) 



5, 



(2) 



-...]) 

S(2) = 



y 

xa,(n-+i(x)n-+i(a;')) + h.c.}. (ElO) 
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Terms which do not have the necessary e^'^^^y structure 
are omitted here. Now we integrate out the fast fields in 
each chain. Then, in chain y we obtain the combination 
like this 

r •or/i</ ^ , /!</ /M /32 r> h l+cos[fe- (03-33^ )] ^ 



(27r)2 



(Ell) 



where X = (x -\- x')/2 and p = x — x' are the center-of- 
mass and relative coordinates. The most relevant term 
that emerges has the spatial derivative acting on the c- 
function which is produced by fast modes (the superscript 
on the integral indicates that it is over the support of the 
fast modes only). In the second line above, the derivative 
has been carried out, bringing down the integral factor 
shown from the derivative of the exponential. After doing 
so, we have approximated the exponential itself by 1. 
This approximation is exact to first order in which is 
infinitesimally small. Then 



Note that 7cone here is itsel f a f unction of run ning R G 
scale ^, as specified in Eq. (E6). Solving Eq. (E16) is 
easy and leads to 



IGttAv 



-4A£ 



-8A£ 



(E17) 



which shows that behavior is controlled by th e 

driving term Jcone/^ right-hand-side of Eq. (E16). 

As discussed above, the chains enter the strongly cou- 
pled limit at C = exp[r] - {v/ where 
Eq. (E17) must be stopped. At this point, the phases 



Oy may be regarded as no longer fluctuating, and hence 
reduce to the classical phases i^y of the main text. Thus 
7bq(-^'0 corresponds directl y to^ bg defined in E q. ( [83] ). 
Combining therefore Eq. (E17) with Eq. (E4), which 



tells us that spontaneous moment of the a-b planes 
^ we arrive at the estimate Eq. (84), ^bq ^ 
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kix sin[fci • p] k2x sin[fc2 • p] 



fk2 



kl 



^2 



X j d2x2cos[2/?(^<-^<+i)]. 



(E12) 



The integral over the relative distance p produces differ- 
ence of two delta-functions, 5{ki^k2) — 5{ki — k2)^ which, 
thanks to the kixk2x factor in the numerator, only dou- 
bles the final result. The whole of the momentum-shell 
integration reduces to 



> ^2 



d k\ 



1 f^' dk dl 



(27r)2 k\ 47r 



47r 



(E13) 



The generated biquadratic correction to the action is 
— 5'(2). Taking 0^ 0^ we see that we indeed generate a 
biquadratic interaction of the form 

i^bq = 7bq ^ dx cos[2P{0y{x) - Oy^i{x))i (E14) 



with 



<^7bq 



cone 



(E15) 



Here we have added a factor v to the generated inter- 
action, accounting for the transformation back to imagi- 
nary time r from xq. The scaling dimension of this term 
is 2 X (2/3)^/(47r) = 8A±. Hence, the RG flow equation 
for 7bq is 



daha = -8A±7bq ■ 



2;;,2 
cone 
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In Sec. |IIC , it was stated that three of the five allowed 
DM couplings, Da^ D'^ and D'^^ can be safely neglected. 
In this appendix, we explain why this is the case. 



Da 



First consider the Da term. As with all the DM cou- 
plings, this is only effective for fields parallel to its D- 
vector, in this case the a axis. For such fields, like the 
Dc term studied in Sec. [VTl it can be "gauged away" for 
a single chain, by an x — dependent spin rotation about 
the z axis of spin. Unlike the Dc coupling, however, the 
Da interaction is constant within each triangular plane. 
Therefore this rotation has negligible effects upon the 
other in-plane couplings, most importantly D'^ = D, 
which we have argued dominates the physics in this field 
orientation, but also J', which plays a subsidiary but 
still important role. This gauge rotation does affect the 
J" interaction, however, since the Da term is staggered 
along the 2: axis, see Eq. (16). We have seen already 



in Sec. IV C that J" itself is already (without Da) in- 



effective in establishing interlayer correlations, and its 
only effects arise through generating the J2 interaction 
between second neighbor layers, Eq. (69). The J2 inter- 
action is, happily, also unaffected by the gauge rotation, 
as the second neighbor layers involved rotate identically. 
Thus even if some Da is present, the analysis of Sec. |IV| 
remains unchanged. 



Dl 



D'a 

interaction, this coupling has the same sign on both di- 



Next consider the D'^ interaction. Unlike the D 



(-^^^) agonal bonds between chains (-D+^ = in Eq. (17)). 
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This means that, hke the J' couphng, this interaction is 
highly frustrated. As a consequence, the leading order 
contributions arising from this term involve a gradient, 
analogous to the twist/cone term in, e.g. Eq. (26). Thus 
the effects of this term are generally strongly suppressed, 
both by this gradient (and associated increased scaling 
dimension) and by its small magnitude, which is of at 
most a few percent. In other words, it carries the same 
scaling dimension as the twist/cone coupling, but is prob- 
ably at least a factor of 10 smaller in magnitude. Thus 
it is always negligible. 



3. D', 

Finally, we turn to the D'^ term. This interaction is 
similar in some ways to the D = D'^ interaction, which 
dominates for fields along a. Both are unfrustrated, as 
they have opposite signs on the two diagonals, and both 
are staggered along the a {z) direction. However, D'^ is 
also staggered along c (?/), while D was constant within 
the triangular planes. 

The fate of D'^ is less clear than that of the prior two 
terms under consideration. It is neither trivially gauged 
away nor obviously negligible. However, it is easy to es- 
tablish that it does compete with many of the key inter- 
actions that have already been identified as driving forces 
in the system. As such, provided D'^ is not too large, it 
loses this competition and has minimal effects. 

First, we see that D'^ has the same scaling dimension 
and hence relevance as the J" term. Moreover, like the D 
term, it competes with the J" interaction because of the 
staggering along z. Hence, if D'^ is not comparable to J" ^ 
it will lose this competition. Indeed, if one assumes the 
form, Eq. (75), which satisfies the 7^ coupling (ex J")^ 
the D'^ term identically vanishes. 

Second, the D'^ term also competes with the Dc term, 
since the latter favors opposite rotations on neighboring 
chains, which the D'^ term attempts to couple. Trans- 
forming to the rotating frame favored by Dc^ Eq. (88), 
will make the D'^ term oscillate, and hence average out 
over long distances. 



Thus to have any significant effect, the D'^ term would 
need to be large enough to overcome at least two com- 
peting interactions. At least for small D'^^ we conclude 
that the phase diagrams established in the main text are 
unchanged. Evidently, this is the case in CS2CUCI4. 



high field magnons with such calculations, the micro- 
scopic paramete rs of th e standard Hamiltonian Eq. (|3| 
were determined) ^^ * ^^ * ^ However, in the standard model, 
all the possible DM vectors are not included. We present 
here a complete analysis based on Eq. (15), which was 
derived in Appendix [B] 

Let us first show that the components of the DM vector 
perpendicular to an applied field can be negligible in the 
spin- wave analysis. To show this, we decompose Si into 
{S) + 5Si^ where {S) is the ordered moment parallel to 
the applied field h. In the linear spinwave theory we 
neglect 5 Si parallel to which means 5 Si x 5Sj is always 
parallel to h and does not couple to the component of 
the DM vector perpendicular to h. One can also show 
that the DM interaction does not produce single magnon 
terms which is proportional to 5Si using the symmetry 
argument in Appendix [B] In what follows, we only retain 
Dq = D ( where ( = h/h. We now take the direction of 
the field {() as a quantization axis of spins. Introducing 
5'f = Si-0 and = ± iS"^ such that = i fj, ( form 
an orthonormal basis, the local Hamiltonian for the bond 
ij is written as 



Hij — JijSi • Sj + Dij^^{S^Sj - S]^Sj) 

= ^(e^^-5+^-+h.c.) + J,,5f5^^, (Gl) 



where Jij = yJfj + D'^j ^ and tan^^j = Dij^^/ Jij. In the 
following, we focus on ^ = ^, ^, and c cases and introduce 



= + (^C)'. tan(/)c = DJJ, (G2) 



j; = J{J'Y tan ^' = D'/J' (GS) 



for on-chain and diagonal bonds, respectively. We now 
apply the Holstein-Primakoff transformation: 

5f = 5 - nu St = {2S - ni)hi, 5" = bl{2S - ni)i , 

(G4) 

with rii = blbi and 5* = 1/2, and obtain 
Hij - JijS{e''*'^^blbj^h.c.)-JijS{ni^nj)^JS^. (G5) 



Appendix G: Spin-wave analysis in a high field 

In this appendix, we study the effect of the DM inter- 
action and the inter-layer interaction J" on the high field 
magnons, and particularly on the ordering wavevector in- 
finitesimally below the saturation field. In a strong mag- 
netic field, the ground state is a fully polarized ferromag- 
netic state and one can easily solve the single-magnon 
problem exactly. By comparing measurements of the 



Denoting by ba,k the Fourier transform of the boson at 
the position R-\-6a^ the spinwave Hamiltonian is written 
as 

Hsw = ^ [H{k) + - (2 J + 4f + 2r)S]^k. (G6) 

k 

where = (^i^^, 62,^, ^3,^, ^4,fc)^ and the 4x4 matrix 
1-L{k) depend on the field direction. 
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1. Field along a axis 

Let us first consider th e cas e of field along the a axis. 
In this case, 1-L{k) in Eq. (|G6|) is given by 



2 \ B, 



(G7) 



Here, the matrices Aa^ki^l^a^ (l^a) ^fc 



/^a,k{^a,(t>'a) 



2Ja COs{kb - (t)a) Jafa{(t>'a'^ k) 



and Bfc = J^\l + e*^'^)/, where / is the 2 x 2 identity 



matrix and fa{(t)'a]k) = (e ^^e^^^ 



0(1 



). We 



also note that /c^ is defined by /c • /i for /i = a, 6, and 
c. We now try to find the location of the minimum of 
the spectrum of 1-magnon excitations, which is given in 
the form of k* = (0,27r(l/2 + e),0). Using the relations 
such as cos(/)a = J/Ja, the eigenvalues of H(0,/c5,0) are 
explicitly obtained as 



^i,± 



J cos ki) + 2 J' cos{k})/2) 



±^{J"Y + {Da sinkb + 2D^^ sin(fe5/2))2, (G8) 
^2,± = J COS kb — 2 J' cos(/c5/2) 

±^{J"Y + [Da Sink, - 2D'^ sin(fc5/2))2. (09) 

Putting = 0, one finds the results consistent with 
Ref. 16, Among the four solutions, uji^± and a;2,±, ^i,- 
has the lowest energy around /C5 = 27r(l/2 + eo), where 
eo = jy (27r J) which is the incommensuration in the ab- 
sence of Da and D'^. It is convenient to introduce a vari- 
able X = cos{ki)/2). Then we rewrite cji^_ as 



2X' 



1 



4- 



4(1 -X2) ( -^X 



J 



(GIO) 



Assuming X is small, we expand the above equation and 
have the following approximate expression for the incom- 
mensuration: 



sin(7re) = 



+ 4(Di)2 



(Gil) 



JV(J")'+4(D;)2 J[{J"Y + 4(D;)2]3/2 



Whether the incommensuraton is enhanced or not de- 
pends on the signs of Da and D'^ and the subtle balance 
of them. We note that, if Da is present, and not too 
small, the incommensuration can be substantially mod- 
ified from the DM- free value, since the second term in 



the brackets in the first line of Eq. (Gil ) is small only in 



the ratio Da/ J' . Since Da was neglected in the experi- 
mental fits in Ref,fT6l, this might lead to small errors in 



the magnetic parameters, at perhaps a level of ten per- 
cent of their estimated values, i.e. an uncertainty in J' 
of iO-1 ^estimated' similarly for D = D'^. Errors of 
the order of 10% of the largest interaction, J, are clearly 
ruled out by the fits.-^ 



2. Field along h axis 

Next we consider the case of field along the h axis. In 
this case, we have 

where 

''^^M={ (G13) 

with h{(l)[]k) = e^^J>(l + e'^^){l + e'^^). We minimize 
the excitation energy to find the ordering wavevector of 
the form /c* = (0,27r(l/2 + e),0). The lowest eigenvalue 
of H(0,/C5,0) is 

uo — J COS kh 



{J"Y - WJ" cos{kt/2) + 4( J^)2 cos2(fc5/2).(G14) 



Similarly to the previous subsection, we rewrite the above 
as 



J = - 1 - ^i:^ - 4:^^ ^ 4 (^^.^ (G15) 

where (J^)^ = ( J')^ + (I^5)^. One can find the incommen- 
suration sin(7re) = —X from the minimum of the above 
equation and observe that DM interaction, D'^^ always 
enhances the incommensuration from its D'^ = value 
J'/{2J). 



3. Field along c axis 

Finally, we consider the cas e of field along the c axis. 
In this case, 'H(fc) in Eq. (G6) is given by 




(G16) 



with fe{4>'e;k) 



2J,cos{kb - ^e) J'JM'c^k) \ 
J'Jc{-<P'c\ -k) 2 Je cos(A:fc + 0^) / 

(G17) 



The minimum of the spectrum of 1-magnon excitation 
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is also of the form A^* = (0,27r(l/2 + e),0). The lowest 
eigenvalue of H(0, k^^ 0) is explicitly obtained as 



uji = -J" ^ J cos h-^{DcY sin^ h + ^{J'Y cos^{h/2). 

(G18) 

Here we have used the relations such as cos 0c = J/Jc- 
The remarkable point here is that the minimum and 
hence the incommensuration e is independent of D^^ and 
J" . So once we know e and J, J', it uniquely determine 
the strength of Dc. Let us now assume that e is of the 
order J^J, which is true if Dc = 0, and obtain approxi- 
mate eigenenergy as 



J 



2sin2(7re) - 2 



sin(7re) 



r 
J 



-1, 

(G19) 



where we have neglected a term proportional to sin^(7re). 
From the above equation, we can obtain the incommen- 
suration e as a function of J, J' and as 



sin(7re) 



2J 



(G20) 



From this relation, we see that the Dc on the J bonds 
enhances the incommensuration e. This is in contrast 
to the measured incommensurability,^^ which is reduced 
compared to the expected one from the ideal standard 
model. 



^ R Anderson, Mater. Res. Bull 8, 153 (1973). 
2 D. A. Huse and V. Elser, Phys. Rev. Lett. 60, 2531 (1988). 
^ B. Bernu, C. Lhuillier, and L. Pierre, Phys. Rev. Lett. 69, 
2590 (1992). 

^ L. Capriotti, A. E. Trumper, and S. Sorella, Phys. Rev. 

Lett. 82, 3899 (1999). 
^ R. Coldea, D. A. Tennant, A. M. Tsvelik, and Z. Tylczyn- 

ski, Phys. Rev. Lett. 86, 1335 (2001). 
^ R. Coldea, D. A. Tennant, and Z. Tylczynski, Phys. Rev. 

B 68, 134424 (2003). 

T. Radu, H. Wilhelm, V. Yushankhai, D. Kovrizhin, 

R. Coldea, Z. Tylczynski, T. Liihmann, and F. Steglich, 

Phys. Rev. Lett. 95, 127202 (2005). 
^ Y. Tokiwa, T. Radu, R. Coldea, H. Wilhelm, Z. Tylczynski, 

and F. Steglich, Phys. Rev. B 73, 134414 (2006). 
^ C. H. Chung, K. Voelker, and Y. B. Kim, Phys. Rev. B 

68, 094412 (2003). 
^° S. V. Isakov, T. Senthil, and Y. B. Kim, Phys. Rev. B 72, 

174417 (2005). 

J. Alicea, O. L Motrunich, and M. P. A. Fisher, Phys. Rev. 
Lett. 95, 247203 (2005). 

M. Kohno, L. Balents, and O. Starykh, J. Phys.: Conf. 
Ser. 145, 012062 (2009). 

M. Kohno, O. A. Starykh, and L. Balents, Nat. Phys. 3, 
790 (2007). 

O. A. Starykh and L. Balents, Phys. Rev. Lett. 98, 077205 
(2007). 

M. Kohno, Phys. Rev. Lett. 103, 197203 (2009). 
R. Coldea, D. A. Tennant, K. Habicht, Habicht, 
P. Smeibidl, C. Wo Iters, and Z. Tylczynski, Phys. Rev. 
Lett. 88, 137203 (2002). 

A. A. Nersesyan, A. O. Gogolin, and F. H. L. Essler, Phys. 
Rev. Lett. 81, 910 (1998). 

M. Bocquet, F. H. L. Essler, A. M. Tsvelik, and A. O. 

Gogolin, Phys. Rev. B 64, 094425 (2001). 

M. Q. Weng, D. N. Sheng, Z. Y. Weng, and R. J. Bursill, 

Phys. Rev. B 74, 012407 (2006). 
2° S. Yunoki and S. Sorella, Phys. Rev. B 74, 014408 (2006). 
2^ Y. Hayashi and M. Ogata, J. Phys. Soc. Jpn. 76, 053705 

(2007). 

2^ T. Pardini and R. R. P. Singh, Phys. Rev. B 77, 214433 



(2008). 

22 H. C. Jiang, M. Q. Weng, Z. Y. Weng, D. N. Sheng, and 
L. Balents, Phys. Rev. B 79, 020409(R) (2009). 

'^^ D. Heidarian, S. Sorella, and F. Becca, Phys. Rev. B 80, 
012404 (2009). 

2^ T. Tay and O. Motrunich, Phys. Rev. B 81, 165116 (2010). 
2^ R. B. Griffiths, Phys. Rev. 133, A768 (1964). 

L Affieck and M. Oshikawa, Phys. Rev. B 60, 1038 (1999). 
2^ N. M. Bogoliubov, A. G. Izergin, and V. E. Korepin, Nucl. 

Phys. B 275, 687 (1986). 

S. Qin, M. Fabrizio, L. Yu, M. Oshikawa, and L Affieck, 
Phys. Rev. B 56, 9766 (1997). 

D. C. Cabra, A. Honecker, and P. Pujol, Phys. Rev. B 58, 
6241 (1998). 

T. Hikihara and A. Furusaki, Phys. Rev. B 69, 064427 

(2004) . 

M. Y. Veillette, J. T. Chalker, and R. Coldea, Phys. Rev. 
B 71, 214426 (2005). 

H. J. Schulz, Phys. Rev. Lett. 77, 2790 (1996). 

A. Kolezhuk and T. Vekua, Phys. Rev. B 72, 094424 

(2005) . 

2^ T. Nikuni and H. Shiba, J. Phys. Soc. of Jpn. 64, 3471 
(1995). 

M. Y. Veillette and J. T. Chalker, Physical Review B 74, 
52402 (2006). 

^'^ E. M. Stoudenmire and L. Balents, Phys. Rev. B 77, 
174414 (2008). 

M. Oshikawa, M. Yamanaka, and L Affieck, Phys. Rev. 
Lett. 78, 1984 (1997). 

A. V. Chubukov and D. L Golosov, J. Phys.: Condens. 
Matter 3, 69 (1991). 

K. Hida and L Affieck, J. Phys. Soc. Jpn. 74, 1849 (2005). 
P. M. Chaikin and T. C. Lubensky, Principles of Con- 
densed Matter Physics (Cambridge University Press, Cam- 
bridge, 1995). 

K. Okunishi and T. Tonegawa, J. Phys. Soc. Jpn. 72, 479 
(2003). 

T. Hikihara, T. Momoi, A. Furusaki, and H. Kawamura, 
Phys. Rev. B 81, 224433 (2010). 

S. Lee, R. K. Kaul, and L. Balents, Arxiv preprint 
arXiv.org:0911.0038 (2009). 



43 



E. F. Shender, Sov. Phys. JETP 56, 178 (1982). 

C. L. Henley, J. Appl. Phys. 61, 3962 (1987). 

M. Yoshida, T. Suzuki, T. Waki, and M. Takigawa, Poster 
presentation at Highly Frustrated Magnetism 2008 confer- 
ence (2008). 

M. Takigawa, private communication (2010). 
M.-A. Vachon, G. Koutroulakis, V. F. Mitrovic, A. Reyes, 
P. Kuhns, R. Coldea, and Z. Tylczynski, J. Phys.: Con- 
dens. Matter 20, 295255 (2008). 
^° T. Ono, H. Tanaka, T. Nakagomi, O. Kolomiyets, H. Mi- 
tamura, F. Ishikawa, T. Goto, K. Nakajima, A. Oosawa, 
Y. Koike, et al., J. Phys. Soc. Jpn. (Suppl.) 74, 135 (2005). 
Z. Weihong, R. H. McKenzie, and R. R. P. Singh, Phys. 
Rev. B 59, 14367 (1999). 

K. Foyevtsova, Y. Zhang, H. O. Jeschke, and R. Valenti, 
J. Phys.: Conf. Ser. 145, 012038 (2009). 
T. Ono, H. Tanaka, O. Kolomiyets, H. Mitamura, T. Goto, 
K. Nakajima, A. Oosawa, Y. Koike, K. Kakurai, J. Klenke, 
et al., J. Phys.: Condens. Matter 16, S773 (2004). 
Y. Fujii, H. Hashimoto, Y. Yasuda, H. Kikuchi, M. Chiba, 
S. Matsubara, and M. Takigawa, J. Phys.: Condens. Mat- 
ter 19, 145237 (2007). 

H. Tsujii, C. R. Rotundu, T. Ono, H. Tanaka, B. Andraka, 
K. Ingersent, and Y. Takano, Phys. Rev. B 76, 060406 
(2007). 

N. A. Fortune, S. T. Hannahs, Y. Yoshida, T. E. Sherhne, 
T. Ono, H. Tanaka, and Y. Takano, Phys. Rev. Lett. 102, 
257201 (2009). 

^'^ J. Ahcea, A. V. Chubukov, and O. A. Starykh, Phys. Rev. 
Lett. 102, 137201 (2009). 
M. Bocquet, Phys. Rev. B 65, 184415 (2002). 
M. Y. Veillette, A. J. A. James, and F. H. L. Essler, Phys. 
Rev. B 72, 134429 (2005). 

D. Dalidovich, R. Sknepnek, A. J. Berlinsky, J. Zhang, and 
C. Kallin, Phys. Rev. B 73, 184403 (2006). 

J. Merino, R. McKenzie, J. Marston, and C. Chung, J. 
Phys.: Condens. Matter 11, 2965 (1999). 
A. E. Trumper, Phys. Rev. B 60, 2987 (1999). 
J. O. Fjaerestad, W. Zheng, R. R. P. Singh, R. H. McKen- 
zie, and R. Coldea, Phys. Rev. B 75, 174447 (2007). 



W. Zheng, J. O. Fjaerestad, R. R. P. Singh, R. H. McKen- 
zie, and R. Coldea, Phys. Rev. B 74, 224420 (2006). 
W. Zheng, J. O. Fjaerestad, R. R. P. Singh, R. H. McKen- 
zie, and R. Coldea, Phys. Rev. Lett. 96, 057201 (2006). 
W. Zheng, R. R. P. Singh, R. H. McKenzie, and R. Coldea, 
Phys. Rev. B 71, 134422 (2005). 
^'^ A. L. Chernyshev and M. E. Zhitomirsky, Phys. Rev. B 
79, 144416 (2009). 

J. Alicea, O. I. Motrunich, and M. P. A. Fisher, Phys. Rev. 
B 73, 174430 (2006). 

J. Alicea, O. L Motrunich, M. Hermele, and M. P. A. 

Fisher, Phys. Rev. B 72, 064407 (2005). 

S. Yunoki and S. Sorella, Phys. Rev. Lett. 92, 157003 

(2004). 

C. Chung, J. Marston, and R. McKenzie, J. Phys.: Con- 
dens. Matter 13, 5159 (2001). 
^2 R. F. Bishop, P. H. Y. Li, D. J. J. Farnell, and C. E. 
CampbeU, Phys. Rev. B 79, 174405 (2009). 
We would like to thank R. Coldea for pointing this exper- 
imental finding to us. 

M. Kenzelmann, R. Coldea, D. A. Tennant, D. Visser, 
M. Hofmann, P. Smeibidl, and Z. Tylczynski, Phys. Rev. 
B 65, 144432 (2002). 
'^^ E. H. L. Essler, A. Furusaki, and T. Hikihara, Phys. Rev. 
B 68, 064410 (2003). 

S. Bailleul, D. Svoronos, P. Porcher, and A. Tomas, C. R. 
Acad. Sci. Ser. 2 313, 1149 (1991). 
T. Moriya, Phys. Rev. 120, 91 (1960). 
'^^ S. Lukyanov and A. Zamolodchikov, Nucl. Phys. B 493, 
571 (1997). 

T. Giamarchi, Quantum physics in one dimension (Oxford 
University Press, 2004). 

I. S. Gradshteyn and I. M. Ryzhik, Tables of integrals, 
Series, and Products (Academic, New York, 6th ed., 2000). 
H. J. Schulz, Phys. Rev. B 34, 6372 (1986). 
O. A. Starykh, A. Furusaki, and L. Balents, Phys. Rev. B 
72, 094416 (2005). 

H. J. Schulz and C. Bourbonnais, Phys. Rev. B 27, 5856 
(1983). 



